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We consider a discrete elliptic equation on the d-dimensional lat- 
tice Z d with random coefficients A of the simplest type: they are 
identically distributed and independent from edge to edge. On scales 
large w.r.t. the lattice spacing (i.e., unity), the solution operator is 
known to behave like the solution operator of a (continuous) ellip- 
tic equation with constant deterministic coefficients. This symmetric 
"homogenized" matrix Ahom = dhom Id is characterized by £ ■ Ab om £ — 
((£ + V0) • A(£ + V</>)) for any direction £ £ R d , where the random 
field (j> (the "corrector" ) is the unique solution of —V* • A(£ + V(j>) = 
such that 4>(0) = 0, Vcj> is stationary and (Vcj>) = 0, (•) denoting the 
ensemble average (or expectation). 

It is known ("by ergodicity") that the above ensemble average of 
the energy density £ = (£ + V<^>) • A(^ + V<f>), which is a stationary 
random field, can be recovered by a system average. We quantify this 
by proving that the variance of a spatial average of £ on length scales 
L satisfies the optimal estimate, that is, var[VJ£77L] < L , where 
the averaging function [i.e., = 1, supp(r/L) C < L}] has to 

be smooth in the sense that |V?7l| < L~ 1 ~ d . In two space dimensions 
(i.e., d — 2), there is a logarithmic correction. This estimate is optimal 
since it shows that smooth averages of the energy density £ decay in 
L as if £ would be independent from edge to edge (which it is not for 
d> 1). 

This result is of practical significance, since it allows to estimate 
the dominant error when numerically computing ahom- 

1. Introduction. 

1.1. Motivation, informal statement and optimality of the result. We 
study discrete elliptic equations. More precisely, we consider real functions 
u of the sites x in a d-dimensional Cartesian lattice Z d . Every edge e of the 
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lattice is endowed with a "conductivity" a(e) > 0. This defines a discrete 
elliptic differential operator —V* • AV via 

-V* ■ (AVu)(x) := £ o(e)(«(x)-u(y)), 

y6^,|x-y|=l 

where the sum is over the 2d sites y which are connected by an edge e = [x,y] 
to the site x. It is sometimes more convenient to think in terms of the 
associated Dirichlet form, that is, 

(Vu • AVu)(x) := v(x)(-V* • (iVu)(s)) 
= - « (2/)) a (e)C"(x) - 

e 

where the last sum is over all edges e, and (x,y) denotes the two sites 
connected by e, that is, e = [x,y] = [y,x] (with the convention that an edge 
is not oriented). We assume the conductivities a to be uniformly elliptic in 
the sense of 

a < a(e) < (3 for all edges e 

for some fixed constants < a < (3 < oo. 

We are interested in random coefficients. To fix ideas, we consider the 
simplest situation possible: 

{a(e)} e are independently and identically distributed (i.i.d.). 

Hence, the statistics are described by a distribution on the finite interval 
[a,/3]. We would like to see this discrete elliptic operator with random co- 
efficients as a good model problem for continuum elliptic operators with 
random coefficients of correlation length unity. 

The first results in stochastic homogenization of linear elliptic equations 
in the continuous setting are due to Kozlov [11] and Papanicolaou and 
Varadhan [18], essentially using compensated compactness. The adaptation 
of these results to discrete elliptic equations in quite more general situations 
than the one considered above (i.e., under general ergodic assumptions) is 
due to Kiinnemann [13] following the approach by Papanicolaou and Varad- 
han for the continuous case, and also to Kozlov [12] (where more general 
discrete elliptic operators are considered). Note that the discrete elliptic op- 
erator —V* • AV is the infinitesimal generator of a random walk in a random 
environment, whence the rephrasing of the homogenization result in [13] as 
the diffusion limit for reversible jump processes in Z d with random bond 
conductivities. With the same point of view, it is also worth mentioning the 
seminal paper by Kipnis and Varadhan [9] using central limit theorems for 
martingales. 
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The general homogenization result proved in these articles states that 
there exist homogeneous and deterministic coefficients Ah om such that the 
solution operator of the continuum differential operator —V • Ah om V de- 
scribes the large scale behavior of the solution operator of the discrete dif- 
ferential operator —V* • j4V. As a by product of this homogenization result, 
one obtains a characterization of the homogenized coefficients ^4hom : it is 
shown that for every direction £ € M. d , there exists a unique scalar field 
(j) such that is stationary [stationarity means that the fields V </>(•) and 
V0(- + z) have the same statistics for all shifts z S and (V</>) = 0, solving 
the equation 

(1.1) - V*-(^(^ + V^)) = inZ d , 

and normalized by (f)(0) = 0. As in periodic homogenization, the function 
7j d 3 x \— > £ • x + 4>{x) can be seen as the yl-harmonic function which macro- 
scopically behaves as the affine function Z rf B x i— > £ • x. With this "corrector" 
0, the homogenized coefficients Ah om (which in general form a symmet- 
ric matrix and for our simple statistics in fact a multiple of the identity: 

^hom — O'hom 

Id) can be characterized as follows: 

(1.2) C-A hom ^=((C + V^-A(C + V4>)). 

Since the scalar field (£ + V0) • A(£ + V</>) is stationary, it does not matter 
(in terms of the distribution) at which site x it is evaluated in the formula 

(1.2) , so that we suppress the argument x in our notation. 

The representation (1.2) is of no immediate practical use, since the equa- 
tion (1.1) has to be solved: 

• for every realization of the coefficients {a(e)} e and 

• in the whole space 7h d . 

In order to overcome the first difficulty, it is natural to appeal to ergodicity 
(in the sense that ensemble averages are equal to system averages), which 
suggests to replace (1.2) by 

(1.3) ^ hom £^]T(£ + V0)-A(£ + V0>7L, 

where rji is a suitable averaging function of length scale L> 1, that is, 

(1.4) supp(7? L )c{|x| <L}, \v L \<L~ d , £>L = 1. 

In fact, on expects the energy density (£ + V<p) ■ A(£ + V0), which is a sta- 
tionary random field, to display a decay of correlations over large distances, 
so that (1.3) seems a good approximation for L S> 1. 

However, one still has to solve (1.1) on the whole space Z rf , albeit for a 
single realization of the coefficients. In order to overcome this second diffi- 
culty, we start with the following observation: since (f> on the ball {\x\ < L} is 
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expected to be little correlated to <\> outside the ball {\x\ > R} for R — L> 1, 
it seems natural to replace <f> in (1.3) by <j)R\ 

(1.5) + V0) • A(£ + V<f>)r)L - + y<t> R ) ■ A{i + Vcf> R ) VL , 

where 4>r is the solution of an equation on a domain (say, a ball) of size R 
with homogeneous boundary conditions (say, Dirichlet): 

- V* • (,4(£ + V0r)) = in Z d n{\x\ <R}, 

(1-6) 

R = O in Z d n{|x| > i?}, 

so that the right-hand side of (1.5) is indeed computable. 

However, V<j) R defined by (1.6) is not statistically stationary, which is a 
handicap for the error analysis. It is therefore common in the analysis of the 
error from spatial cut-off to introduce an intermediate step which consists 
in replacing equation (1.1) by 

(1.7) T- 1 (j} T - V* ■ {A{£ + V(f> T )) = inZ* 

Clearly, the zero order term in (1.7) introduces a characteristic length scale 
y/T (the notation T that alludes to time is used because T" 1 corresponds 
to the death rate in the random walker interpretation of the operator T _1 — 
V* • AV). In a second step, (1.7) is then replaced by 

T^Vr - V* • (A(£ + V0r,fl)) =0 in Z d n {\x\ < R}, 

(j>T,R = in Z d n{|x| > R}. 

The Green's function Gt(x,]j) of the operator T _1 — V* • AV is known to 
decay faster than any power in \^ y \ <C 1 uniformly in the realization of the 
coefficients [see, in particular, Lemma 2.8(iii)]. Therefore, one expects that 
4>T and 4>t,r agree on the ball {\x\ < L} up to an error which is of infinite 

order in e = jfrz ( e ^ s * ne i nverse of the distance of the ball {\x\ < L} to 
the Dirichlet boundary {\x\ = R} measured in units of vT, see, e.g., [2], 
Section 3, for related arguments). Hence, we shall consider + V4>t) • 
A(£ + V4>t)vl as a very good proxy to the practically computable + 

V(t>T,R)-A(C + V<P T ,R)VL- 

Y,(C + V0r) • A(£ + V0 T )riL « + V<fr,R) • ^(C + ^^t,r)vl- 

In view of this remark, we restrict our attention to the error we make 
when replacing 

It is natural to measure this error in terms of the expected value of its 
square. This error splits into two parts, the first arising from the finiteness 
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of the averaging length scale L and the other arising from the finiteness of 
the cut-off length scale VT: 



(1.2) 



+ V0 T ) • A(£ + V<Ar)r? L - ((£ + V0) • A(£ + V$> 



(1.8) = var[^(£ + V^ T )-A(£ + Vc/> T )7 ?L 

(1.9) + | (]T(£ + V</>r) • A(£ + V0 T )r/ L ) - <(£ + V0) • A(£ + V$> 

In view of the stationarity of (f + V0 T ) ■ A(£ + V0r), of (1.4) and of (1.1), 
the second part (1.9) of the error can be rewritten as 

+ V</> T ) • A(£ + V<fr)-q L ) - <(£ + V0) • A(£ + V0)> |* 

(i.io) = i((e + v^t) • A(e + v<at) - (e + v</>) • + V0)>| 2 

= ((V(j)T - V(/>) • A(V0T - V</>)) 2 . 

What scaling can we expect for the two error terms (1.8) and (1.10)? A 
heuristic prediction can be easily inferred from the regime of small ellipticity 
contrast, that is, 1 — % <C 1 (and a = 1 w.l.o.g.). In this regime, to leading 
order, the two error terms (1.8) and (1.10) behave like 



var 



and (|V0r-V0| 



2\2 



x\52(t-(A-(A))Z + 2Z.V$)TiL 
where 4> and 4>t are defined via 

(1.11) -A$ = V*-((A-(A))Q, 

(1.12) T- 1 T -A0T = V*-((A-(A))O, 

respectively. In the first error term, we have replaced (j>T by </> for simplicity 
of the exposition. 

These error terms can be computed in a straightforward manner. Indeed, 
as shown in the Appendix, they scale for any direction |£| = 1 as: 

(1.13) var[]T(£-(yl-<yl))£ + 2£-V0)7 ?L ] ~ L~ d 

( T~ d , for d < 4, 

(1.14) {\V<fr- V^| 2 ) 2 ~ I T- 4 ln 2 T, ford = 4, 

{T~ 4 , ford>4. 

We now argue that the first error term (1.13) is the dominant one (in 
dimensions d < 8). In order to do so, we argue that the choice of L ~ y/T 
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is natural [for which (1.13) dominates (1.14) in dimensions d< 8]. Indeed, 
we recall that in the ball {\x\ < L], 4>t is a proxy for the computable 4>t,r 
(defined on the larger ball {\x\ < R}). The error is of infinite order in the 
distance between the two balls, measured in the length scale VT, that is, in 
e := vT/ (R — L) <C 1. Hence, for the sake of discussing rates, we may indeed 
think of L ~ VT ~ R. 

In this paper, we therefore focus on the error term (1.8) coming from 
the finite range L of the spatial average. In Theorem 2.1 (see also Remark 
2.1), we shall establish that (1.13) holds as an estimate also for its nonlinear 
counterpart (1.8), that is, 

(1.15) var[^(£ + V0 T )-A(£ + V0 T )r?L 

with two minor restrictions: 



• In dimension d = 2, the prefactor depends logarithmically on T (whereas 
for dy^2, the prefactor depends only on the ellipticity constants). 

• The spatial averaging function tjl has to be smooth in the sense that 
\Vr] L \ <L- d ~ l in addition to (1.4). 

The estimate for the higher order term (1.9) will be the object of a subse- 
quent work. 

1.2. Discussion of the works of Yurinskii and of Naddaf and Spencer. In 
this subsection, we comment on two papers on error estimates (in the sense 
of the previous subsection) which from our perspective are the essential ones. 
We also explain how our work relates to these two papers. 

Still unsurpassed is the first quantitative paper, the inspiring 1986 work 
by Yurinskii [21]. He essentially deals with the error (1.9) arising from the 
spatial cut-off T. In our discrete setting of i.i.d. coefficients a(e) and for 
dimension d > 2, his result translates into 

(1.16) (\V<fr - V0| 2 ) < T (2-d)/(4+d)+,5 ) 

for T>1 and some arbitrarily small 5 > 0, see [21], Theorem 2.1 (and [5], 
Lemma A. 5, for this rephrasing of Yurinskii's result). 

Yurinskii derives estimate (1.16) by fairly elementary arguments from the 
following crucial variance estimate of the spatial averages ^tVl of §t on 
length scales L: 

(1.17) varpT^J <T\-j) 

for 1< T < L d and some arbitrarily small 5 > 0, see [21], Lemma 2.4. Let 
us comment a bit on the proof of (1.17): by stationarity of 4>t, the variance 
can be reformulated as a covariance, that is, 



var |^ (j) T rj L = cov ^ <f) T fj L ; <j) T (0) 
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with a modified averaging function fj^. The starting point for (1.17) is to 
control the covariance by: 

(i) An additive decomposition of 4>t(0) over all finite subsets S of the lattice 
Z d , that is, 4>t{0) = ^2scz d <fcr,s(0), where 4>t,s(0) only depends on aig, 
that is, the coefficients a restricted to the subset S. 

(ii) An estimate on how sensitively Y14>tVl depends on aig. 

The decomposition in (i) is based on the probability measure on path space 
[0, oo) 3 t 1 y r](t) £ Z d describing the random walk generated by the operator 
— V* • AV (for a fixed realization of a). Indeed, this probability measure on 
path space allows for a well-known representation of c/>t(0) in terms of paths 
starting in (via the expected value). Hence, the splitting can be obtained 
from restricting the expected value to all paths n with image S (up to some 
exit time larger than T), see [21], Lemma 2.3. 

The sensitivity estimate (ii) comes in form of the deterministic energy- 
type estimate 

|5Z0t??l-^^l| 2 <^ (l + |V^ T (e)| 2 ), 

edges e s.t. eC\Sy^0 

where 4>t is the solution of T -1 ^ — V* ■ A(£ + V4>t) = with coefficients 
A which differ from A only on the subset S, see [21], (1.17). 

The third ingredient for (1.17) is an estimate of the probability that a path 
n starting in crosses a given edge e. This probability can be estimated in 
terms of the Green's function Gt{x, 0) of the operator T _1 — V* • AV (where 
x is one of the two sites on the edge e). Yurinskii then appeals to estimates 
on Gt(x,u) that only depend on the ellipticity bounds a < a < j3 of A (and 
therefore do not depend on the realization of a) see [21], Lemma 2.1. As is 
well known, these type of estimates rely on the Harnack inequality. 

Our variance estimate (1-15) also relies on these deterministic estimates 
of the Green's function Gt(x,v), see Lemma 2.8. However, our strategy to 
estimate a variance differs substantially from Yurinskii's strategy of (i) and 
(ii). As a matter of fact, with our methods, we could derive the optimal 
variance estimate 



(1.18) vax\52<hriL 



2-d 



for L>1. Estimate (1.18) is optimal in the sense that we obtain the above 
scaling in the regime of "vanishing ellipticity ratio" 1 — ^ ^ 1 by the argu- 
ments in the previous subsection. Still, the optimal estimate (1-18) would 
not yield the optimal estimate (1-14) by Yurinskii's argument to pass from 
(1.17) to (1.16). 

Our strategy of estimating a variance is inspired by an unpublished paper 
by Naddaf and Spencer [17]. They use a spectral gap estimate to control the 
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variance of some function X of the coefficients {a(e)} e dges e (i-e., a random 
variable) : 

(1.19) var[X]</ V 



/ dX 



v edges e 



\da(e) 



see [17], page 4. This type of estimate can be seen as a Poincare estimate 
with mean value zero w.r.t. the infinite product measure that describes the 
distribution of the coefficients (and the optimal constant in this estimate 
is given by the smallest nonzero eigenvalue of the corresponding elliptic 
operator, whence "spectral gap"). Naddaf and Spencer derive (1.19) via 
the Brascamp-Lieb inequality for a large class of statistics for {a(e)} e d ge s ei 
which however does not include all i.i.d. statistics of {a(e)} e dges e considered 
by us. We therefore rely on a slight modification of (1.19), see Lemma 2.3. 

We also follow Naddaf and Spencer in the sense that we treat the variance 
of an energy density. However, they express their result not in terms of the 
energy density of 4>t but of a generic solution u with a compactly supported, 
deterministic right-hand side /, that is, 

(1.20) -V* • AVu = V* ■ f. 

Using (1.20), they obtain the formula ^ u ' AVu = — | Vu(e)| 2 so that 

an application of (1.19) yields the following estimate on the energy density 
X = £ Vu ■ AVu: 



[1.21) varl^Vu-AVu <(^2\Vu\ 



see [17], Proposition 1. 

Naddaf and Spencer also remark that provided the ellipticity contrast 
1 — % is small enough, Meyer's estimate holds which states that 

(i-22) Ei Vu i 4 ^Ei/i 4 ' 

with a constant that only depends on a, (5. The combination of (1-21) and 

(1.22) yields the a priori estimate 

(1.23) var[EVn-AVn] <]T|/| 4 > 

see [17], Theorem 1. Since the left-hand side of (1.23) scales as (volume) 2 , 
while the right-hand side only scales as volume, this estimate reveals the op- 
timal decay of fluctuations on the macroscopic level, very much like (1.15). — 
There is a somewhat theatrical convention in the homogenization literature 
to call the lattice spacing e instead of 1 which highlights this scaling. Fol- 
lowing Naddaf and Spencer, we use Meyer's estimate, albeit applied on the 
Green's function Gx{x,y), see Lemma 2.9. 
We will make use of the following notation: 
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• d > 2 is the dimension; 

• J" zd dx denotes the sum over x S Z d , and f D dx denotes the sum over 
x £ Z d such that x £ D, D open subset of R d ; 

• (•) is the ensemble average, or equivalently the expectation in the under- 
lying probability space; 

• var[-] is the variance associated with the ensemble average; 

• < and > stand for < and > up to a multiplicative constant which only 
depends on the dimension d and the constants a, (3 (see Definition 2.1 
below) if not otherwise stated; 

• when both < and > hold, we simply write ~; 

• we use 3> instead of > when the multiplicative constant is (much) larger 
than 1; 

• (ei, . . . , e<i) denotes the canonical basis of 
2. Main results. 

2.1. General framework. 

Definition 2.1. We say that a : 7L d x 7L d — > M + , (x,y) h> a{x, y) is a con- 
ductivity function on Z d if there exist < a < f3 < oo such that: 

• a(x,y) =0 if \x-y\ ^ 1, 

• a(x, y) = a(y, x) E [a, /3] if |x - y\ = 1. 

We denote by ^4q,^ the set of such conductivity functions. 

Definition 2.2. The elliptic operator L:Lf oc (Z d ) -)• Lf oc (L d ) 
associated with a conductivity function a 6 „4. Q/ g is defined for all x 6 Z rf by 

(2.1) (Lu)(z) = -V*-A(x)Vu(x), 

where 



Vu(x) : = 

and 



it(x + ei) — u(x) 
it(x + ej) — u{x) 



X7*u(x) :-- 



u{x) — u(x — ei) 
u(x) — u{x — e^) 



A(x) := diag[a(x, x + ei), . . . , a(x, x + e^)]. 
In particular, it holds that 

(Lu)(x)= ^2 a ( x iy)( u ( x )- u (y))- 

y,\x-y\=l 

If a(x,y) = 1 for \x — y\ = 1, then the associated elliptic operator L is the 
discrete Laplace operator, and is denoted by — A. 
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Definition 2.3 (Discrete integration by parts). Let d > 2, he L 2 (Z d ) 
and g G L 2 (Z d ,M. d ). Then the discrete integration by parts reads 

/ h(x) V* • g(x) dx = — / Vh(x) ■ g(x) dx. 
Ji d Jz d 

We now turn to the definition of the statistics of the conductivity function. 

Definition 2.4. A conductivity function is said to be independent and 
identically distributed (i.i.d.) if the coefficients a(x,y) for \x — y\ = 1 are 
i.i.d. random variables. 

Definition 2.5. The conductivity matrix A is obviously stationary in 
the sense that for all z G Z d , A{- + z) and A(-) have the same statistics; and 
for all x, z G Z d , 

(A(s + z)) = (A(x)). 

Therefore, any translation invariant function of A, such as the modified 
corrector 4>t (see Lemma 2.2), is jointly stationary with A. In particular, 
not only are 4>t and its gradient V0t stationary, but also any function of 
A, 4>t and V0t- A useful such example is the energy density (£ + V(/>t) • 
A(£ + V0t)j which is stationary by joint stationarity of A and V(pT- 

Another translation invariant function of A is the Green functions Gt 
of Definition 2.7. In this case, stationarity means that Gt{- + z,- + z) has 
the same statistics as Gt(-,-) for all z G so that in particular, for all 
x,y,z G Z d , 

(G T (x + z,y + z)) = {G T (x,y)). 

Lemma 2.1 (Corrector ([13], Theorem 3)). Let a G A a p be an i.i.d. con- 
ductivity function, then for all £ G M. d , there exists a unique random function 
cj) : Z rf — > M which satisfies the corrector equation 

(2.2) - V*- A(x){V(f>(x)+£) = inZ d , 

and such that 0(0) = 0, Vcj) is stationary and (Vcft) = 0. In addition, (\ V</>| 2 ) < 

We also define an "approximation" of the corrector as follows. 

Lemma 2.2 (Approximate corrector ([13], proof of Theorem 3)). Let 
a G A a p be an i.i.d. conductivity function, then for all T > and £ G W 1 , 
there exists a unique stationary random function 0t : Z rf — )• 1R which satisfies 
the "approximate" corrector equation 

(2.3) T~ 1 (l) T {x)-V t ■ A(x){V(h{x) + C) = U inZ d , 
and such that (4>t) = 0. In addition, T -1 + (|V0 T | 2 ) < |£| 2 - 
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Note that <f>x is stationary, whereas is not. 

Definition 2.6 (Homogenized coefficients). Let a G A a p be an i.i.d. 
conductivity function and let ^ £ W 1 and </> be as in Lemma 2.1. We define 
the homogenized d x d-matrix j4hom as 

(2.4) £ • A h om£ = <(£ + V0) • A(£ + V$(0)>. 

Note that (2.4) fully characterizes A-^ ora since -Ahom is a symmetric matrix 
(it is in particular of the form ahomld for an i.i.d. conductivity function). 

2.2. Statement of the main result. Our main result shows that the energy 
density £ := T" 1 ^ + (V^t + • A(V4>t + °f the approximate corrector 
4>t, which is a stationary scalar field, decorrelates sufficiently rapidly so that 
smooth spatial averages (defined with help of fluctuate as they would if 
£ would be independent from site to site (as is the case for the tensor field 
A of the coefficients). The strength of fluctuation is expressed in terms of 
the variance. In more than two space dimensions (i.e., d > 2), the estimate 
does not depend on the cut-off scale \/T and thus carries over to the energy 
density of the corrector (p. In two space dimensions, we are not able to rule 
out a weak (i.e., logarithmic) dependence on the cut-off scale \/T: 

Theorem 2.1. Let a 6 A a/ 3 be an i.i.d. conductivity function, and let <p 
and 4>t denote the corrector and approximate correctors associated with the 
conductivity function a and direction £ G M. d , |£| = 1. We then define for all 
L > and T ^> 1 the symmetric matrix A^^ characterized by 

i-A L ^:= [ (T-VT(x) 2 + (V0 T (x)+e)-^)(V0 T (x)+e))r/L(x)dx, 
Jz d 

where x i— > t]l(x) is an averaging function on (—L,L) d such that J^ d r/i(x) dx = 
1 and ||V?7l||£cx) <L~ d ~ 1 . Then, there exists an exponent q>0 depending 
only on a, (3 such that 

ford = 2 w^[i-A L ^]<L- 2 {\nTy, 

(2.5) 

for d > 2 var[£ • A L rf] < L~ d . 

In particular, for d>2, the variance estimate (2.5) holds for the energy 
density of the corrector (ft itself. 

Remark 2.1. While it is natural to include the zero-order term T~ 1 (^) 
into the definition of the energy density, it is not essential for our result. Here 
comes the reason: by a simplified version of the string of arguments which 
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lead to Theorem 2.1 we can show that the variance of the zero-order term 
is estimated as 



< / QnT)«, for d = 2, 
L 2 ~ d , for d > 2. 



Hence, this term is of lower order in the regime (of interest) L < T. 

The main ingredient to the proof of Theorem 2.1 is of independent inter- 
est. It states that all finite stochastic moments of the approximate corrector 
(j>T are bounded independently of T for d > 2 and grow at most logarithmi- 
cally in T for d = 2. 

Proposition 2.1. Xei a G «4 a/ g be an i.i.d. conductivity function, £ G K d 
wii/i |£| = 1 and Zei </>t denote the approximate corrector associated with 
the conductivity function a, and £. Then there exists a continuous function 
7 : M + — > R + suc/i traat for all q G M + , there exists a constant C q such that 
for all T > 0, 

ford = 2 (\M0)\ q )<C q (lnT)^\ 

(2.6) 

ford>2 (\M0)\ q )<C q . 
In addition, j(2n) = n(n + 1) for all n = 2 l , I G N large enough. 

Let us give a heuristic argument for the behavior of (I^t^O)! 9 ) for d=l. 
In this case, for T = oo, the gradient of the corrector associated with £ = 1 
is explicitly given by 

V0=-t4k-1. 

a(a 1 ) 

Hence, ^(x) G M behaves as a discrete Brownian motion in x G Z once we 
have fixed its value at 0. Usually, one imposes 0(0) = almost surely, so 
that for |x| ~ VT, 

M*)H<Ko)=o)~(v^r /2 . 

Yet, one may choose a nontrivial initial value. In particular, one may also 
consider </>(0) = 0t(O) (which yields a corrector field different from the one 
in Definition 2.1). With cj) defined this way, 4>t{x) and <f>(x) are expected not 
to differ much provided \x\ <C \fT. On the one hand, from this we deduce 
that <f>r(x) behaves locally as a discrete Brownian motion starting at </>t(0), 
so that we have as above 

<|<M*)-<Mo)H~M 9/2 
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for all q > and \x\ <C VT. On the other hand, since §t is stationary, 

(\M*) - <M0)H < <|<M*)I 9 > + (|^(0)| 9 ) = 2<|<M0)H. 
These two estimates indeed suggest that 

<i<Mo)n>v^/ 2 -, 

where the minus sign accounts for the fact that the argument only holds 
for \x\ <C VT — we may for instance miss logarithmic corrections. Hence, 
there is a transition between unboundedness and boundedness in T for some 
d£ (1,3). The linearization of the problem in the regime of vanishing ellip- 
ticity contrast, that is, 1 — % -C 1, suggests that d = 2 is indeed the critical 
dimension for Proposition 2.1, that is, the dimension where a logarithmic 
behavior is to be expected. However, there is no reason why d = 2 should 
be critical for Theorem 2.1. Indeed, in the case of d = 1, the statement of 
Theorem 2.1 holds without a logarithm. 

In view of our discussion of the case d = 1 and the observations in case 
of vanishing ellipticity contrast, it is not surprising that the statement of 
bounded stochastic moments is harder to prove the closer we are to d = 2. 
For the experts in homogenization, let us give a quick sketch of the strategy 
of the proof of this result. Independent of the dimension, the proof always 
starts from the variance estimate (Lemma 2.3) applied to ^(O) 9 and makes 
use of the representation of 9 g^^ with help of the gradient V x Gt(x,0) 
(Lemma 2.4). 

• In the case of d > 4, the uniform pointwise, but suboptimal, decay | V x Gt(x, 
y)\ < \x — y\ d ~ 2 , which can be easily obtained from the same pointwise de- 
cay of the Green's function itself, is sufficient. 

• In case d = 4, it would be enough to appeal to the Holder estimate (with 
exponent 7 only depending on the ellipticity contrast) in order to get the 
somewhat better pointwise decay \V x GT{x,y)\ < \x — y\ d ~ 2 ~"' ■ 

• In d = 3, we need (in addition) the optimal decay \V x Gt(x, y)\ < \x — 

which cannot be a pointwise control, but only an average control 
on dyadic annuli. In fact, we need the control of the square average, which 
we easily obtain from the Cacciopoli estimate. 

• For d = 2, the square average is not sufficient anymore, we need the average 
to some power p > 2, as provided by Meyers' estimate (Lemma 2.9). This 
forces us — somewhat counterintuitively — to first estimate high moments 
of 4>t, so that the exponent we put on the gradient of the Green's function 
can be chosen close to 2 (and thus below Meyers' exponent). 

In this presentation, we only display the last strategy (although it is an 
overkill for dimensions d > 2). 

As a corollary of Proposition 2.1, we obtain the following existence and 
uniqueness result of stationary solutions to the corrector equation (1.1) for 
d > 2, which settles a long-standing open question. 
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Corollary 2.1. Leta£A a /3 be an i.i.d. conductivity function. Then, 
for d>2 and for all £ S M. d , there exists a unique stationary random field <j> 
such that ((f)) = and 

-V* • (A(£ + V0)) = inZ d . 
In addition, {(f) 2 + \V<p\ 2 } < |£| 2 . 

We will not prove Corollary 2.1 in detail. Here comes the argument. 
Proposition 2.1 yields the a priori estimate {4>j) < C which is uniform in 
T. This additional estimate allows us to pass to the limit in the probabil- 
ity space for <fir, as it is done for V<^r m [13], proof of Theorem 3. Note 
that the corrector fields of Lemma 2.1 and Theorem 2.1 do not coincide 
(only their gradients coincide). Uniqueness further requires the argument 
by Papanicolaou and Varadhan in [18], which does not appear in [13]. 

Let us point out that Proposition 2.1, Theorem 2.1 and Corollary 2.1 
hold true for more general distributions, provided the variance estimate of 
Lemma 2.3 below holds. In particular, the law of a(x, x + e^) may depend on 
the direction e^, which would give a general diagonal homogenized matrix 
(not necessarily a multiple of the identity matrix). More generally, a(x,x') 
and a(y, y') may also be slightly correlated. We do not pursue this direction 
in this article. 



2.3. Structure of the proof and statement of the auxiliary results. Not 
surprisingly, in order to control the variance of some function X of the coef- 
ficients a (like the spatial average of the energy density of the approximate 
corrector </>t), one needs to control the gradient of X w.r.t. a. As in [17], 
this is quantified by the following general variance estimate: 



Lemma 2.3 (Variance estimate). Let a = {oj}j g N be a sequence of i.i.d. 
random variables with range [a, (3] . Let X be a Borel measurable function of 
a £ R N (i.e., measurable w.r.t. the smallest a-algebra on M N for which all co- 
ordinate functions K N 3 a i— > aj G M. are Borel measurable, cf. [10], Definition 

114). 

Then we have 

2\ 



(2.7) 



var[X] < ( Y^sup 



dX 



da; 



var[ai], 



where sup a . \ ^-\ denotes the supremum of the modulus of the ith partial 
derivative 



ox 

da; 



of X with respect to the variable a% S [a, f3] . 
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Remark 2.2. Let us comment a bit on Lemma 2.3. Estimate (2.7) is a 
weakened version of a spectral gap estimate 

(2.8) 

\ 1=1 

which already played a central role in Naddaf and Spencer's analysis of 
stochastic homogenization [17], Section 2. We note that for i.i.d. random 
variables, such a spectral gap estimate (2.8) follows "by tensorization" from 
the one-dimensional spectral gap estimate 

(2-9) {X{a^)-{X{ ai )f<(\^Xy 

see, for instance, [14], Lemma 1.1. The one-dimensional spectral gap estimate 
(2.9) holds under mild assumptions on the distribution of a\. Yet, (2.9) does 
not hold for atomic measures like (X(ai)) = ^(X(l) +X(2)). Since Lemma 
2.3 covers the case of atomic measures, we only obtain the weaker form (2.7) 
of (2.8). Despite this technical detail, the proof of Lemma 2.3 is very similar 
to the one in [14], Lemma 1.1. 

As in [17], in the proof of Theorem 2.1, we will make use of the fact 
that T~ l (f>T + {S7(j)T + " A(\7<pT + is an energy density, which yields the 
following elementary formula for the partial derivative w.r.t. the value a(e) 
of the coefficient in the edge e = [z, z + e^]: 

(T-Vf + (V£r + ■ A(V<p T + 0)(x)vl(x) dx 



da(e) 

(2.10) = -2 J (J^^Vl • A{V<t> T + £)) (x) dx 

up to minor modifications coming from the discrete Leibniz rule, see Step 1 
of the proof of Theorem 2.1. 

This formula makes the gradient of the averaging function r]L appear; in 
order to benefit from this, we assume that the averaging function is smooth 
so that we get an extra power of The merit of (2.10) is that we need 
to control the partial derivative 9 q^^ of the approximate corrector 4>t( x ) 
(and not of its spatial derivatives). Not surprisingly, this partial derivative 
involves the Green's function Gt(x,-). More precisely, it involves the gra- 
dient V Zi GT(x, z) of the Green's function with singularity in z [and not its 
second gradient V Zj S7 x Gt{xi z), for which we would not have the optimal 
decay rate uniformly in a]. 

We define discrete Green's functions as follows. 
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Definition 2.7 (Discrete Green's function). Let d > 2. For all T > 0, 
the Green's function Gt '■ A a p x Z d x Z d — > Z d , (a, rc, y) i— )• Gt(x, y; a) associ- 
ated with the conductivity function tt is defined for all y G and a G 
as the unique solution in L^(Z rf ) to 



(2.11 



T 1 GT(x,y;a)v(x) dx 



+ f Vv(x)-A(x)V x G T (x,y;a)dx = v(y) Vv£L 2 (Z d ), 
where A is as in (2.1). 

Note that the existence and uniqueness of discrete Green's functions is a 
consequence of Riesz' representation theorem. Throughout this paper, when 
no confusion occurs, we use the short-hand notation Gt(x, y) for Gt(x, y; a). 

The following lemma provides the elementary formula relating the "sus- 
ceptibility" 9 Qaff °f ^t( x ) t° the Green's function GT(x,y). 

Lemma 2.4. Let a G A a p be an i.i.d. conductivity function, and let Gt 
and (fir be the associated Green's function and approximate corrector for 
T > and £ G M d , |£| = 1. Then, for all e=[z,z + e;] and x G Z d , 

(2.12) dMx,a) = _ {ViMz . a) + Gt ( z x . o)) 

oa{e) 

and for all n G N, 



sup 

a(e) 



9 -[<AT(x;«r +1 ] 



<9a(e) 

(2.13) <\Mx-,a)\ n (\ViMz-,a)\ + l)\V Zi G T (z,x-a)\ 

+ (|V^ T (z;a)| + l) n+1 |V Zi G T (z,x;a)r +1 . 
In addition, it holds that 

(2.14) sup|Vi£r(z;a)| < |V^ T (z;a)| + 1. 

a(e) 

iVoie £/ia£ t/ie multiplicative constant in (2.13) depends on n next to a, (3 
and d. 

In addition, Lemma 2.4 provides uniform estimates on 8 ^(e)* " * n a ( e ) 
(the case n > 1 is needed in Proposition 2.1). In order to obtain this uniform 
control in a(e), we need to control V 2 G(z,x;a) uniformly in a(e). Again, 
this comes from considering ^'gj^p^ ■ The following lemma provides the 

elementary formula for ^^q}^^ an d a uniform estimate in a(e). 
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Lemma 2.5. Let Gt ■ A a/ 3 xZ^xZ^l, (a, x, y) \-t Gt{x, y, a) be the 
Green's function associated with the conductivity function a for T > 0. For 
all e = [z,z + ei\ and for all x, y G Z d ; it ZioWs £/ia£ 

(2.15) -^-jG T {x, y; a) = -V Zi G T (x, z; a)V Zz G T (z, y; a). 

As a by-product, we also have: for all x G Z rf 

(2.16) sup|V 2i G r (z,x;a)| < |V^Gr(z,ac;o)|. 

o(e) 

There is a technical difficulty arising from the fact that a has infinitely 
many components. In Lemma 2.3, this technical difficulty is handled by the 
strong measurability assumptions on X. The following lemma establishes 
these measurability properties for (fxr, so that we can apply Lemma 2.3. 

Lemma 2.6. Let a G A a p be an i.i.d. conductivity function, and let 
Gt(-,-;cl) and (Pt{']o) be the associated Green's function and approximate 
corrector for £ G M. d , d>2, and T > 0. Then for fixed x, y G Z rf , Gt(x, y, •) 
and <Pt(x;-) are continuous w.r.t. the product topology of A a p (i.e., the 
smallest/ coarsest topology on M. E , where E denotes the set of edges, such 
that the coordinate functions W E B a \— > a e G M. are continuous for all edges 
e£E). 

In particular, Gy(x,y;-) and 4>t(x;-) are Borel measurable functions of 
a G A a /3, so that one may apply Lemma 2.3 to •) and nonlinear funtions 
thereof. 

The proof of Theorem 2.1 crucially relies on the fact that 4>t is almost 
bounded independently of T (in d > 2). More precisely, it relies on the 
fact that any moment (0y(O) n ) is bounded independently of T as stated 
in Proposition 2.1. Starting point for Proposition 2.1 is again Lemma 2.3, 
which is iteratively applied to <^r(0) m where m increases dyadically. This 
is how Lemma 2.4 comes in again. However, the crucial gain in stochastic 
integrability is provided by the following lemma. It can be interpreted as a 
Cacciopoli estimate in probability and relies on the stationarity of 4>t- 

Lemma 2.7. Let a G A a fi be an i.i.d. conductivity function, and let 4>t 
be the approximate corrector associated with the coefficients a for £ G M. d , 
|f | = 1. Then for d>2 and for all n G 2N, 

(2.17) <|<Mo)r(|v<Mo)| 2 + |v*<Mo)| 2 )> < <l^r(o)>, 

where the multiplicative constant does depend on n next to a, (3, and d, but 
not on T > 0. 



18 



A. GLORIA AND F. OTTO 



In order to prove Proposition 2.1 via Lemma 2.3 [applied to </>;r(0) n ] and 
Lemma 2.4, we need some weak version of the optimal decay of the gradient 
V z Gt(x, z) of Green's function in \x — z\, that is, 

(2.18) \\7 z Gt{x, z; a)\ < \x — ~ d uniformly in a and T. 

This decay is the best we can hope as can be checked on the Green function 
for the Laplace equation. The same decay property is needed to prove The- 
orem 2.1 via Lemma 2.3 [applied to (2.10)] and Lemma 2.4. Yet it is well 
known from the continuum case that there are no pointwise in z bounds 
of the type (2.18) which would hold uniformly in the ellipticity constants 
a, f3. (An elementary argument shows that any bound on V x G(x,y) which 
would be uniform in a and inl/2<|x — y|<l would yield that a bounded 
a-harmonic function has bounded gradient. However, for d = 2 and for any 
7 > 0, there are examples of a-harmonic functions from the theory of quasi- 
conformal mappings that are not Holder continuous with exponent 7, see 
[6], Section 12.1.) Nevertheless, (2.18) holds in the square averaged sense on 
dyadic annuli, as can be seen by a standard Cacciopoli argument based on 
the optimal decay of the Green's function itself, that is, 



in the case d > 2. The pointwise estimate (2.19) in x and z is a classical result 
[7], Theorem 1.1, that relies on Harnack's inequality. It has been partially 
extended to discrete settings, see in particular the Harnack inequality on 
graphs [3]. However, we did not find a suitable reference for the BMO-type 
estimate in the case of d = 2. On the other hand, we do not require the 
pointwise version of (2.19), but just an averaged version on dyadic annuli. 
The statements we need are collected in the following lemma. 

Lemma 2.8. Let a S A a p, T > and Gt be the associated Green's func- 
tion. For all d>2 and q> 1, r > 0, 

(i) BMO and L q estimate: for all i?» 1, 



(2.19) 



Gt(x,z) <\x- z 



2-d 



uniformly in a and T, 



(2.20) ford = 2 




q dx<R\ 



(2.21) ford>2 




where GT(-,y){\ y - x \<R} denotes the average o/Gr(-,y) over the ball 
{x £ Z d , \x-y\< R}. 
(ii) Behavior for R ~ and d = 2: 
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(iii) Decay at infinity: for all R > \[T , 

(2.23) f G T {x,y) q dx<R\R 2 - d ) q {VfR- 1 ) r . 

J R<\x-y\<2R 

The multiplicative constants in (2.20), (2.21) and (2.23) depend on q,r next 
to a, (3 and d. 

We present a proof of Lemma 2.8 which for d = 2 is a direct version of 
the indirect argument developed in [4], Lemma 2.5, in case of a nonlinear, 
continuum equation. For the convenience of the reader, we also include the 
proof for d > 2 — anyway, it has the same building blocks as the argument 
for d = 2. This makes our paper self-contained w.r.t. the properties of Gt- 

However, it is not quite enough to know (2.18) in the sguare-averaged 
sense on dyadic annuli. In order to compensate for the fact that we only 
control finite stochastic moments of V<^(0) via Proposition 2.1, we need to 
control a pth power of the gradient V z Gt{x,z) of Green's function in the 
optimal way for some p > 2. This slight increase in integrability is provided 
by Meyers' estimate, which yields such a p > 2 as a function of the ellipticity 
bounds a, j3 only. Meyers' estimate has already been crucially used in [17], 
however in a somewhat different spirit. There it is used that for sufficiently 
small ellipticity contrast, 1-I<1, one has p > 4. The following lemma is 
the version of Meyers' estimate we need and will prove. 

Lemma 2.9 (Higher integrability of gradients). Let a £ A. a p be a conduc- 
tivity function, and Gt be its associated Green's function. Then, for d>2, 
there exists p > 2 depending only on a, (3, and d such that for all T > 0, 
p>q>2, k>0 and i? » 1, 

(2.24) f \V z G T (z,0)\ q dz<R d (R 1 - d ) 9 mm{l,VTR- 1 } k . 

J R<\z\<2R 

For technical reasons, we need a pointwise decay of GT(x,y;a) in \x — y\ 
uniformly in a (but not in T). The decay we obtain is suboptimal and easily 
follows from Lemmas 2.8 and 2.9 using the discreteness. 

Corollary 2.2. For all d>2 and T > 0, there exists a bounded radially 
symmetric function hx £ L x (% d% ) depending only on d,a,(3, and T such that 

G T {x,y;a) <h T (x-y) 

for all x, y £ 7L d and a E A a p. 

Lemmas 2.8 and 2.9 only treat Gt away from the diagonal x = y — which 
is a consequence of the fact that the scaling symmetry is broken by the 
discreteness. Using the discreteness, the following corollary establishes a 
bound independent of T and a. 
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Corollary 2.3. For all a e A a /s, T > and x,y £ 

\VG T (x,y;a)\<l. 

Finally, for the proof of Theorem 2.1, we need to know that also the 
convolution of the gradients of the Green's functions decays at the optimal 
rate, that is, 

I \V z G T (x,z)\\V z G T {x\z)\dz 
Jz d 

(2-25) 

< \x — x | uniformly in a and T. 

As for (2.18), it is not necessary to know (2.25) pointwise in (x,x'), but 
only in an averaged sense on dyadic annuli. The following lemma shows that 
(2.25) for linear averages can be inferred from (2.18) for quadratic averages. 



Lemma 2.10. Let h T E Lf oc (Z d ) be such that for all 1 and T > 0, 

(2.26) ford = 2 / h%{z) dz < min{l, VTR' 1 } 2 , 

JR<\z\<2R 

(2.27) ford>2 [ h 2 T (z) dz < R 2 ~ d , 

J R<\z\<2R 



and for R ~ 1 

(2.28) ford>2 / h 2 T (z)dz<\. 

J\z\<R 

Then for R^> 1 



(2.29) 



ford = 2 I I hT{z)hT{z — x)dzdx 
'|x|<R. 



<i? 2 max{l,ln(v / T J R -1 )}, 
(2.30) ford>2 / / h T (z)h T (z - x) dz dx < R 2 . 

J\x\<RJZ d 

We present the proof of Proposition 2.1 and Theorem 2.1 in Section 3. We 
gather in Section 4 the proofs of the decay estimates for the discrete Green 
functions (i.e., Lemmas 2.8 and 2.9, and Corollaries 2.2 and 2.3) since they 
are needed at several places in the paper, and may be of independent inter- 
est. The proofs of the remaining auxiliary lemmas are the object of Section 5. 
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3. Proofs of the main results. 

3.1. Proof of Proposition 2.1. Starting point are Lemmas 2.3 and 2.6, 
which yield 

var[<M0n<E( su P 

e \ a(e) 



da(e) 



where ^ e denotes the sum over the edges. Using now (2.13) in Lemma 2.4, 
this inequality turns into 

var[<M0n< / X;(0r(O) 2(m - 1) (|V^ T (^)| + l) 2 |V, i G' T (^O)| 2 

+ (\ViMz)\ + lf m \V Zz G T (z,0)\ 2m ) dz, 

where we have replaced the sum over edges e by the sum over sites z£Z (! 
and directions for i £ {1, . . . , d} according to e = [z, z + e*]. Simplifying 
further, we obtain 

var[0 T (O) m ]< / (M0) 2im - 1} (\VMz)\ + l) 2 \V z G T (z,0)\ 2 

(3.1) 

+ (|V^ T W| + l) 2m |V,G T (z,0)| 2m )^. 

We proceed in four steps. Assuming first that for n big enough and for 
all m < n it holds that 

I {M^ 2{m ~ 1] {\^Mz)\ + if\^ z G T (zM 2 

JZ d 

(3.2) +(\V<p T {z)\ + l) 2m \V z G T {zM 2rn )dz 

< {{M 0)^ r /n-l/Hn + X)) + 1} | InT, for d = 2, 

we prove the claim in the first step. The last three steps are dedicated to 
the proof of (3.2) for n large enough. 

Step 1. Proof that (3.1) and (3.2) imply (2.6). 

For notational convenience, we set fJ*d{T) = 1 for d > 2 and /J-d{T) = InT 
for d = 2. Let n = 2 l , I G N*. Using the elementary fact that 

<<M0) 2m ) < (0T(O) m ) 2 + var[^ T (O) m ], 

from the cascade of inequalities (3.1) and (3.2) for m = 2 l ~ q , q € {0, ...,/}, 
we deduce 

<<Mo) 2 - 2i > < <<Mo) 2 ') 2 + ^((Mof^-vwn+i)) + 1} 
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(estimate 0) 

Mo) 2 ' 2 " 9 ) < <<M0) 2i " 9 ) 2 + Mr)((M0) 2 Y /2? - 1/(re(re+1)) + 1) 

(estimate q) 

(Mo) 2 ' 2 °) < (Mo)) 2 + fi d (T)((Mo) 2n ) 1/n - 1/(n{n+1)) + 1) 



Lemma ^.2^ 



(estimate I). 



We then take the power 2 q of each (estimate q) and obtain using Young's 
inequality: 



(0) 2n ) < <M0) n > 2 + // d (T)(<0 T (O) 2n ) 1_1/(n(n+1)) + 1), 



(Mo) 2 - 2l -Y<(Mor q ) 2q+1 

(3 ' 3) (MorT +1 <(Mor q - 1 ) 2q+2 

+ ^ T f + \(Mo) 2n ) 1 - 2q+1/{n{n+1)) + i), 



(M0f) n < fJ.d(T) n ({(pT(0) 2n ) 1 ~ 1/in+1) + 1). 

Since the multiplicative constants in each line of (3.3) only depend on 
a,/3,d,n and q, a linear combination of these I + 1 inequalities with suit- 
able positive coefficients allows us to cancel the respective terms both on 
the left- and right-hand sides, which yields 

l 

(3.4) (Mo) 2n ) <^2^ d (Tf\(<hm 2n ) 1 ~ 2q/(nin+1)) + 1)- 

Using Young's inequality, each term gives the same contribution and (3.4) 
turns into 

(3.5) <Mo) 2n ><^en n(n+1) . 

Formula (2.6) is then proved for all q < 2n using Holder's inequality in 
probability. 



VARIANCE ESTIMATE FOR EFFECTIVE DIFFUSION 



23 



Step 2. Estimate for the Green's function. 

Let p > 2 be as in Lemma 2.9. We shall prove that for all q > 1 and R S> 1 
the following holds 

ford = 2 f \V z G T (z,0)\ q dz 

JR<\z\<2R 

(3.6) 

<R 2max ^ q ^R- q mm{l,VTR- 1 } q , 
ford>2 / \V z G T (z,0)\ q dz 

JR<\z\<2R 

< ^dmax{l,g/p}/j^l-d\g 



(3.7) 



We split the argument into two parts to treat q > p and q <p, respectively. 
For q > p, we use the discrete L p — L q estimate: 

r \ 1 / q ( f \ 

/ \V z G T (z,0)\ q dz) < / \V z G T (z,0)\ P dz) . 

JR<\z\<2R J \JR<\z\<2R J 

Combined with (2.24) in Lemma 2.9, it proves (3.6) and (3.7). 

For q<p, we simply use Holder's inequality with exponents (|,^r^) in 
the form 

R~ d [ \V z G T (z,0)\ q dz) ^ < (r"* [ \V z G T {zM P dz\ '\ 

JR<\z\<2R J \ JR<\z\<2R J 

that we also combine with (2.24). 
Step 3. General estimate. 

Let % > be a random variable. In order to prove (3.2), we will need to 
estimate terms of the form 



( X \V z G T (z,0)\ q ) 1/r dz 

i 

for q,r > 1. Relying on (3.6) and (3.7), we show that 

(xIv.Gt^h 1 /^ 



(3.8) 

if dmaxi 1, 1 -- + — [> + (1 - d)- < 0, 



< (x) 1/r 



d>2, 



r rp 



InT, if 2max<i 1,1-- + — 1 --=0, 

r rp J r 



d = 2. 

Note that there is no overlap in (3.8). For d > 2, we will only make use of 
the estimate with dmaxjl, 1 — - + ^-} + (1 — d)^ < 0. For d = 2, we will use 
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the estimate both with 2max{l, 1 — - + ^-} — £ < 0, and with 2max{l, 1 — 

^ + ^-} — 2 = 0, which requires a specific argument. 

Let i m i n 6 N,i m ; n ~ 1 be such that Lemma 2.9 holds for R > 2 lmin . To 
prove (3.8), we use a dyadic decomposition of Z d in annuli of radii R L = 2 l : 

( X \V z G T (z,0)\ 9 ) 1/r dz 

i 

(3.9) = f ( X \V z G T (z,0)\ q ) 1/r dz 

J|z|<2 I min 



oo „ 

+ E / ( X |V z G r (z,0)r) 1 /'-dz. 
i-i . JRi<\z\<R i+1 

t> — I'm in 



IR l <\z\<R i+1 

Using Corollary 2.3, we bound the first term of the right-hand side by 

(x\V z G T (z,0)\ q ) l/r dz<(x) 1/r . 



'|*|<2*n 

For the second term of the right-hand side, we appeal to Holder's inequality 
with (r, 



oo „ 

E / ( X \V z G T (z,0)\ q ) 1/r dz 

i=iini JR><\z\<R l+1 



>R l <\z\<R i+1 

1/r 



< E ^ 1/r (l 



{ X \VzG T (zM q )dz 

Ri<\z\<Ri +1 
so that (3.9) turns into 

( X \V z G T (z,0)\ q ) 1/r dz 

\V z G T (z,0)\ q dz ) 

'Ri<\z\<R i+1 



= (x) 1/r + E ( J Rf) 1 " 1/r /x / |V*G T M)|*dz 
. \ JRi<\z\<R i+1 

Using then (3.6) and (3.7), we get 



X I \V z G T (z,0)\Uz 

Ri<\z\<Ri+i 



< 



( X )R^ max{1 ' q/p} R; q min{l, VTRr l y, d = 2, 
(x)R- max{hq/p} (Rt d ) cl , d>2. 



Hence, 

( X \V z G T (z,0)\ q ) 1/r dz 



< 



VARIANCE ESTIMATE FOR EFFECTIVE DIFFUSION 25 

oo 

? 2max{l,l-l/r+g/(rp)}-g/r . r 1 .fifr>-l\q/r 



( X ) 1/r ^maxll.l-l/r+g/^j-g/r ^-1}?/^ d = ^ 

i=0 
oo 

(x) 1/r ^ i2 dmax{1,1 ~ 1/r+9/(r?,)}+(1_d)9/r d > 2 



i=0 



We distinguish two cases. If (imax{l, 1 — i + ^-} + (l — d)^<0, then 

oo 

/ <x|V,G T M)|*> 1/r < ^^HW-iWWWi^ < (x) iA. 

This proves the first estimate of (3.8). For d = 2, and 2max{l, 1 — - + — 
- = 0, then we obtain 

r 1 



/ ( X \V z G T (z, 0)\^/ r dz < (x) 1 ^ Vmin{l, y/TR^} 

/ oo 

<( X ) 1/r lnT + ^i?, 



\ 8=0 / 

<( X ) 1/r (l+lnT). 

This proves the second estimate of (3.8). 
Step 4. Proof of (3.2). 

Let n > 1 and n >m> 1. We first treat the first term of the left-hand 
side of (3.2). In that case Holder's inequality in probability with exponents 
(n + 1, and the stationarity of V^t show 



(^(o^-^dv^^l + iflVA^o)! 2 )^ 

i 

< I (<|v^)l 2(n+1) ) 1/(n+1) + i) 

(3.10) x (|0 r (O)| 2 ^ 1 ^ n+1 )/ n |V,G T (z,O)| 2 ( n+1 )/ n ') ri/(n+1) ^ 

= ((|V<AT(0)| 2(n+1) ) 1/(n+1) + l) 

x / (|^T(0)| 2(m ^ 1)(ri+1)/n |V 2 G r (z,0)| 2 ( ri+1 )/ n ') ri/(n+1) dz. 



We apply Lemma 2.7 to bound the first ensemble average in (3.10): 

<|v<Mo)l 2(n+1) > 



,i=l 
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(3-11) d s 

stationant y /^ |v ^ (0)|2 ^ (0)2 J 
\i=l / 

(2.17) 



We now want to apply Step 3 to the right-hand side integral of (3.10), that 

... 2(n+l) 

is, setting q = y ' 
involves the number 



is, setting q = ^±^, r = 2±i and % = |(M0)| 2(m_1)(n+1)/n - Estimate (3. 



(3.12) 



dmaxj 1, 1 h — \ + (1 - d)- 

[ r rp ) r 

= dmaxi 1, — 1 — i + 2(1 — d). 

[ n + l p) 



We distinguish the cases d > 2 and d = 2. For d > 2, we have that the number 
(3.12) is equal to d + 2(1 — d) = 2 — d and thus negative for n sufficiently 
large since p > 2. Hence, (3.8) yields 

(|</> T (0)| 2(m - 1)(n+1)/n | V z G T (z, o)! 2 ^ 1 )/™)"/^ 1 ) dz 

< (\4 >T (o)\ 2 ( m - 1 )( n + i y n ) n /( n+ v < (^(o)! 2 ™)^ 1 ^, 

where in the last inequality we appealed to Jensen in probability using 

2(m-l)(n + l) < 2(n-l)(n + l) < ^ 
n n 

The combination of this with (3.10) and (3.11) yields as desired 

{Mu) 2{m - l) (\vMz)\ + tfW z G T {zM 2 )dz 

< ^^o) 2n ) 1 /( n+1 )+( ?n - 1 )/ n + 1 = (^(o) 2n ) m / n ~ 1 /( n ( n+1 ^ + 1. 

We turn to the case d = 2. We note that the number (3.12) is zero for n 
large enough since p > 2. Thus, from (3.8), we infer as we did above that 

{Mv? {m ~ l) {WMz)\ + rf\VzG T (zM 2 )dz 

d 

< (lnT)(((/>(0) 2n > m/ri - 1/{n(n+1)) + 1). 

Let us now treat the second term of the left-hand side of (3.2), which differs 
from the first term only when m > 2. As for the first term, Holder's inequality 
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in probability with n "^+1 )> * ne stationarity of V<^t and Lemma 2.7 

imply 



((|V0T(^)| + ir m |V Zi G T ( z ,O)|^}^ 

i 

< f {{\vMz)\ 2{n+1) ) m/{n+1) + i) 

(3.13) x (|V 2 G T (^,0)| 2 ( n+1 ) m /( n - m+1 V n " m+1)/(n+1) ^ 

< «<Mo) 2n > m/(n+1) + 1) 

x f (|V,G T (z,0)| 2 ( n+1 ) m /( n - m+1) ) (n " m+1)/(n+1) ^. 



We use (3.8) with X = 1, 9 = ^m+T an< ^ r = n.-m+i > ^ n wmcn case we have 
(3.14) 



dmax<i 1, 1 -- + —} + (1 - d)~ 



Tfl 2 777, 

dmax<! 1, 1 }> + (1 - d)2m. 

n + 1 p 



We claim that this number is negative for n sufficiently large. Indeed, if 

m I 2m 
n+1 p 



max{l, -™ +22*} = !, then 



dmax 1, 1 } + (1 - d)2m = d + 2m(l - d) 

[ n + 1 p J 

= (2m-l)(l-d) + l<0 
since d > 2 and m > 2. Otherwise, max{l, ^ I + 2 f} = ^ l + 2 f, and 

dmaxjl,— ^— + — 1 + (1 - d)2m = 2m( d( — — — - + - ) + 1 - d 
I n+1 p J V ; V \2{n+l) pj 

< 2m(l - < 

for d > 2 and n large enough since - <\- This shows that (3.14) is negative 
so that we obtain by (3.8) 

{\V Z G T (Z, )\^+l)m/(n-m+l)^n-m+l)/(n+l) ^ < L 

! 

Combining this with (3.13) yields 

((|v<M*)l + i) 2m 1 v,g t (z, o)| 2m ) dz < (Mo) 2n ) m/(n+1) + 1 
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= (^ T (0) 2n ) m / n ~ m//(n ( n+1 )) + 1 



< (</> T (0) 2n ) m//n ~ 1/(n(n+1)) + 1. 
This concludes the proof of the proposition. 

3.2. Proof of Theorem 2.1. Let us define the spatial average of a function 
h : Z rf — > M with the mask t)l by 

((h)) l ■= / h(x)rj L (x)dx, 

where rji satisfies 

5 d ->[0,l] supple) c(-L,L) d , 



(3.15) 



/ r) L (x)dx = l, |Vr7i|<L" 



The claim of the theorem is that there exists q depending only on a, f3, and 
d such that 

var^T" 1 ^ + (V0t + ■ A(V<fr + £)»l] < L- d fj, d (Ty, 

where fJ-d{T) = 1 for d > 2 and Hd(T) = InT for d = 2. Since we are not 
interested in the precise value of q, we adopt the convention that q is a 
nonnegative exponent which only depends on a, (3, and d but which may 
vary from line to line in the proof. 

Starting point is the estimate provided by Lemmas 2.3 and 2.6 



(3.16) 



var[((T- + (V0r + £) ■ ^4(V0 T + 



< 



e a ( e ) 



(9 



((T-^T + + • ^(V0T + 0))l 



da(e) 



Step 1. In this step, using the notation e= [z,z + ej, we establish the 
formula 







da(e) 



(3.17) 



«T- Vr + (V£r + £) • ^4(V0 T + C)))x 
= 2 / (V i 0r^)+ei)V* t Gr(«,x) 

x ^^a(x-e j ,x)Vfa( x )(y*<l >T (x) 
+ r /L (z)(V^ T + ^) 2 ^). 
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Indeed, by definition of we have 

^-{(T-VI + (V^t + 6 • A(V<fr + £)»z 

= / ri L (x)— (T- 1 tf T + (V0 T + O-A(y<p T + dx. 
J% d da(e) 



We note 



da(e) 



so that 



aa(e) aa(e) 

+ (V0t + O-^tt(V0t + O (x) 
aa(e) / 

+ (v^T+ei) 2 w^-z), 



* -((r- VI + (V^r + £M(V£r + £)>>£ 



<9a(e) 

(3 ' i8) - 2 1 (" £ ( T ~ v + v Hi • A(v ^ + «)) (i) ^ 

+??L(2)(v^ T +e i ) 2 ^)- 

Using the discrete integration by parts formula of Definition 2.3, the first 
term of the right-hand side of (3.18) turns into 

(3.19) =-/ Pf-{x)V*-{ri L A{V<h + m*)te 

J Z d da(e) 

+ JA VLT ~^ T ~^) ix)dX - 
We now use the following discrete Leibniz rule: 

V* • (vLAiVfa + 0)(x) = V l(x)(V* ■ A(Vcf>T + 0)(x) 

d 

+ E ^*Mx)[A{V<t> T + ZMx - ej), 
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where [A(VcpT + 0]j denotes the jth coordinate of the vector A(V(pT + £)• 
For notational convenience, we take advantage of the diagonal structure of 
A (although this is not crucial) to rewrite the latter equality in the form 

(3.20) =VL(x)(V*-A(V<j>r + Z))(x) 

d 

i=i 

The combination of (3.20) with (3.19) and the use of the equation satisfied 
by 4> T , 

T~Vt - V* • A(V(j)T + = 0, 

yield 

\l ( T-'ch Ife + V |fe • A(V4>t + 0) ) (*) dx 



d(pT 
d da(e) 



da(e) da(e) 
d 



» - ej,x)V*7] L (x)(V*(l)T(x) + ^j) j dx. 



Using now Lemma 2.4, this turns into 



v^( T "^4) +v ^y A(V0T+e) 



(3.21) 



(2.12) 



(V i «Mz) + &)Vs i G r (*,s) 



x ^^a(x-e j ,x)V*r ?L (x)(V*^ T (^) + Ci)J dx. 

Inserting (3.21) into (3.18) proves (3.17). 

Step 2. In this step, we provide the estimate 

d 



sup 

a(e) 



da(e) 



(3.22) 



< 



«r- Vt + &<fr + ■ A{v<p T + 0))l 



|V 2 G T (z,rE)||V*r ?L (x)|(|V^ T (x)| 2 + |V<M^)| 2 + l)dx 



+ 7 ?L (z)(|V0 T (^)| 2 + l). 
Indeed, from Step 1, the boundedness of a, and |£| = 1, we infer that 
8 



sup 

a(e) 



da(e) 
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(3.23) 



< f (sup|Vi0r(z)| + l)sup|V 2i G r (z,x)||V*77 L (x 

JZ d V a(e) 7 o(e) 



x I sup | V*^r(a;)| + 1 ) 

V a(e) ' 

+ r ?L (z)fsup|V^ T (^| 2 + l). 

V a(e) y 

Hence, in the remainder of this step, we have to deal with the suprema 
over a(e). Recalling that e = [z,z + ei], the two following inequalities are 
consequences of Lemmas 2.5 and 2.4: 

(2.16) 

sup|V Zl G T (z,x)\ < \V Zi G T (z,x)\ forallxeZ d , 

a(e) 

(2.14) 

sup|Vi^r(z)| < |V<0r(z)| + l. 

a(e) 

The last inequality we need is 



(2.14) 



sup|V*^r(x)|<|V*0 T (x)|+sup|V i 0T(z)| + l < |V^ T (x)| + |V^ T (z)| + l. 

a(e) a(e) 

It is then proved combining the boundedness of a and the following bound 
on the derivative of V*0t(^) with respect to a(e): 







da(e 



■V*0t(x) 







da(e 



(2.12) 



|V*((V^ T (^) + e 4 )V, l G T (z,x))| 
\(ViM*) + Si)VziV*xGT(z,x)\ 



< 2(|V^ T (z)l + l&l) ™p |VG T | 

Z d xZ d 

< |V^ T (^)| + 1, 

where we have used the uniform bound on VGt provided by Corollary 2.3. 
Combining these three inequalities with (3.23) yields 

d 



sup 

a(e) 



da(e] 



< 



(\V</>t(z)\+1)\V*Gt(z,x)\\V*til(x)\ 

i 

x (|VVr(s)| + |V0r(*)l + l)<k 
+ r /L (z)(|V^ T (z)| 2 + l) 
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from which we deduce (3.22). 

Step 3. In this step, we argue that 

var[«T-V! + (Vcf)T + ■ A{V<fr + £)})z] 

(3.24) <([ (I |V z G T (z,x)||V*77 L (x)||V*<M*)| 2 ^ dz 

\Jz d \Jz d J 

(3.25) +([([ \VzG T (z,x)\\V* VL (x)\\VMz)\ 2 dx) dz 

\Jz d \Jz d J 

(3.26) +11(1 \V z Gr{z,x)\\V*ri L {x)\da\ dz 

\Jz d \Jz d J 

(3.27) +(/ ri L (zf(\V<h(z)\ 2 + l) 2 dz 

\Jz d 

Indeed, inserting (3.22) in (3.16) yields 
var[«T-VT + (V0t + ■ A(V4> T + £))> L ] 



< 



^{L l V * G H^)l|VV(*)l(|V*Mz)| 2 + |V0 T (z)| 2 + l)dx 

+ (5> 2 >)(|v<Mz)| 2 + i; 



We then use Young's inequality in the first term of the right-hand side of 
this inequality and we replace the sum ^ e over edges [z, z + e»] by d times 
the sum over z £ Z d to establish this step. 

It now remains to estimate the terms (3.24), (3.25), (3.26) and (3.27) to 
conclude the proof of the theorem. 

Step 4. Estimate of (3.27): 

(3.28) U VL {zf(\VMz)? + I) 2 dz\ < fid(T) q L~ d . 

Indeed, by stationarity we have 

d 

<|VM*)| 4 > < + e,)| 4 + \4>t{z)\ A ) = 2d(M0) 4 ), 

i=l 

so that 

VL (z) 2 (NMz)\ 2 + rfdz) <( [ VL (z)\\VMz)\ 4 + 1) dz 

z d 

r, L (z) 2 ((\VMz)f) + l)dz 
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<«<M0) 4 > + 1) / VL(z) 2 dz. 
Jz d 

On the one hand, it follows from Proposition 2.1 that 

<<M0) 4 }<^(T)*, 
with q = 7(4). On the other hand, it follows from (3.15) that 



/ VL (z) 2 dz<L- d . 
Jz d 



This establishes Step 4. 
Step 5. Estimate of (3.26): 

(3.29) n \V g Gr{z,x)\\V*ri L (x)\dx^ dz^ < Md (T)« L~ d . 

We expand the square 

1(1 \V z G T {z,x)\\V* VL (x)\dx) dz\ 

\V*i lL {x)\\V* VL {x')\\V z G T {z,x)\\V z G T {z,x')\dxdx' dz 



VV(x)||VV(x')l / (\V z GT(z,x)\\V z G T (z,x')\)dzdxdx'. 



We then use Cauchy-Schwarz' inequality in probability and the stationarity 
of G T : 

(\V z G T (z,x)\\V z G T (z,x')\) 

<(\V x G T (z,x)\y/ 2 (\V z GT(z,rf)\Y 2 

= (\V z G T (z - x, U)\ 2 ) l l 2 {\V z G T (z - x'M 2 ) 1 ' 2 . 

Hence, with the notation 

h{y) :=(\V y G T (y,0)\ 2 ) 1/2 , 
we have by definition of rj^: 

f (( \V z G T (z,x)\\V* VL (x)\dx) dz) 

< I |V*r ?L (x)||V*r ?L (x / )| / h(z - x)h(z - x')dzdxdx' 



< L -2(d+i) I h{z-x)h{z-x')dzdxdx' 

J\x\<L J\x'\<L . 



34 A. GLORIA AND F. OTTO 

= L -2(d+l) [ [ [ h{z')h{ Z ' +X- X ')dz' dxdx' 

J\x\<L J\x'\<L JZ d 

<L- d ~ 2 [ [ h(z')h(z' -y)dz'dy. 

J\y\<2L JZ d 

We note that 

h 2 (y)dy=( [ \V y G T {yM 2 dy). 

*R<\y\<2R \JR<\y\<2R I 

On the one hand, for ]J>1 we have according to Lemma 2.9 (for q = 2) 
ford = 2 f h 2 (y)dy<R 2 - 2 mm{l,VTR- 1 } 2 

JR<\v\<2R 



lR<\y\<2R 

= min{l,v / T-R~ 1 } 2 , 
f r>drjyl-d\2 

lR<\y\<2R 



ford>2 / h 2 (y)dy<R d (R 1 
JR<\v\<2R 



R 2 ~ d . 



On the other hand, for R ~ 1, Corollary 2.3 implies 

ford>2 / h 2 (y)dy<l. 
J\y\<R 

Hence, we are in position to apply Lemma 2.10, which yields as desired 



\<2L 



h(z')h(z' -y)dz'dy<L 2 Li d (T). 



Note that for d = 2,we have used the elementary fact that max{l, In \[TL 1 } < 
InT. 

Step 6. Estimate of (3.25): 



\V z G T (z,x)\\V* VL (x)\\VMz)\ 2 dx) dz)<n d {TYL- d . 

\JZ d \JZ d 

(3.30) 

As in Step 5, 

7 z G t (z,x)||V*7? l (x)||V0t(2)| 2 ^ dz 



d \Jz d 

V*v L (x)\\V*r] L (x')\ 



d Jz d 



[ (|V M T (z,x)\\V z G T (z,x')\\V(j) T (z)\ 4 ) dz dx dx' . 
Jz d 
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This time, we use Holder's inequality with (p,p, in probability (where 
p > 2 is the exponent in Lemma 2.9): 

(|V z G t (z,x)||V,G t (z,x')||V^t(^)| 4 > 

<(|V,Gt(z,x)| p ) 1 / p (|V,Gt(z,x')| p ) 1/p (|V^ t (^)| 4p/{ ^ 2) ) (p " 2)/p - 
By stationarity of Gt and c^r, we obtain with Proposition 2.1 

(\V z G T (z,x)\\V z G T (z,x')\\VMz)\ 4 ) 

< ^ d {T)"{\V z G T {z - x, 0)\P) 1/p (\V z G T (z - x', 0W) 1/P - 
Hence, with the notation 

h(y) :=(\V y G T (y,0)\ P ) 1/P , 

by definition of r\^\ 

\V z G T (z,x)\\V*r) L (x)\\V4> T (z)\ 2 dx\ dz 



Vd{T) q I I |V*j7l(x)||V*t/l(^)I / h(z-x)h(z-x')dzdxdx' 

z d 



<li d {TYL~ d - 2 [ [ h(z')h(z'-y)dz'dy. 

J\y\<2L JZ d 

As in Step 5, we shall establish that for R 3> 1 

ford = 2 f h 2 (y)dy<mm{l,VfR- 1 Y 
J R<\y\<2R 

(3.31) 

ford>2 / h\y)dy<R 2 - d , 

JR<\y\<2R 

and for R ~ 1 

(3.32) ford>2 [ h 2 {y)dy<\. 

J\y\<R 

Once this is done, Lemma 2.10 implies as desired 



\<2L 



h(z')h(z' - y) dz dy < L Hd(T), 



using in addition that max{l, In y/TL -1 } < InT for d = 2. As above, (3.32) 
is a direct consequence of Corollary 2.3. We now deal with (3.31). Note that 
according to Lemma 2.9, we have for R 3> 1 



(3.33) 



ford = 2 f h p (y)dy<R 2 ~ p mm{l,VfR^ 1 } p , 

JR<\y\<2R 



ford>2 f h p (y)dy<R d (R l - d ) p . 

JR<\v\<2R 



( R- d I h 2 {x)dx\ 

\ JR<\x\<2R J 



< ( R~ d 
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We now argue that this yields (3.31). Indeed, by Jensen's inequality 

v 1/2 

h?(x)da 

>R<\x\<2R 

\ Vp 

hP{x)dx 

R<\x\<2R J 
(3.33) r ( jR -2 fi 2-p min { 1) ^ jR -l}P ) l/P ) d=2] 

1 (A- d i2 d (J? 1 - d )*») 1 /p, d>2, 

Rr 1 min{l, VTR- 1 }, d = 2, 
R^ d , d>2, 

which implies (3.31). 

Step 7. Estimate of (3.24): 

(3.34) n n^ z G T {z,x)\\V*nL{x)\\V*Mx)\ 2 dx\ dz\ < n d (T)« L 
As in Steps 5 and 6, 

\ 2 

|2 



| V z G T (z, x) 1 1 V* VL (x) 1 1 V* 4> T (x) I l dx J dz 

\V*r) L (x)\\V* V L(x')\ 

x f {\V z G T (z,x)\\V z G T (z,x')\ 
Jz d 

x \V*ct) T (x)\ 2 \V*(P T {x')\ 2 }dzdxdx'. 

Holder's inequality with (p,p, ^^r^) ™ probability (where p> 2 is 
exponent in Lemma 2.9) then yields 

(\v z g t (z, x )\\v z g t (z,x')\\v*m^\ 2 \^*M^)\ 2 ) 

<(\V z G t (z,x)\p) 1 /p(\V z G t (z,xT) 1/p 

x (|VVT(x)| 4p /^- 2 ))^ 2)/(2p \|V*0 T (xO| 4p/ ^ 2) ) b_2)/(2p) . 
The stationarity of Gt and (j>T, and Proposition 2.1 show 

(\V z G T (z,x)\\V z G T (z,x')\\V*Mx)\ 2 \V*Mx')\ 2 ) 

< ii d (T)*(\V s G T {z - x,0)\ P ) 1/p (\V z G T (z - x', 0)\n 1/p - 
We may now conclude as in Step 6. 



VARIANCE ESTIMATE FOR EFFECTIVE DIFFUSION 37 

The theorem follows from the combination of Step 3 with (3.28), (3.29), 
(3.30) and (3.34). 

Step 8. Extension to the energy density of the corrector field for d > 2. 
Let AL tOC be defined by 

6 • A L ^ := / (V<Xx) + • A(x)(V<f>(x) + 0^l(x) dx, 
Jz d 

for all L^> 1. The claim is 

var[£-A ii00 £]<L- d , 
for d > 2. It is proved using (2.5) provided we show 
(3.35) var[£ • A LiO0 £] < liminf var[£ • A L>T (\. 

As we shall prove, the following two convergences hold: 

(S + V<fr{x)) ■ A(x)(£ + V<fr{x))ii L (x) dx 



(3.36) (^ + V0(x))-A(x)(^ + V^(x))/i L (x)dx 
\Jz d 

which in fact amounts to ((£ + V<j) T ) • + V^r)) -)• ((£ + V</>) • A(£ + V<£)) 
by stationarity, and 

(£ + V0 T (x)) • A(x)(£ + V0 T (x))^ L (x) dx 

! 

(3.37) / (e + V^(x))-A(x)(C + V^(x))/x L (x)dx 
H d 

weakly in probability. 

We may now conclude the proof of (3.35). Expanding the variance, one has 



var[£ • A L , T £] = / (J (S + VM%)) ■ Mx)(Z + V<fr{x))n L (x) dx 

(£ + V0 r (x)) • A(x)(£ + VMx))^l{x) dx) . 
! / 

By (3.36), the second term of the right-hand side converges to (£ • ^4hom£) 2 
as T — > oo, whereas by lower-semicontinuity of quadratic functionals, (3.37) 
implies that 

(f + V0(x)) • A(x)(£ + V<t>(x))n L {x) dx 

1 j 

< liminf/ (I (f + V<fr{x)) • + V(/. T (x))/i L (x) dx 
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which shows (3.35). 

It remains to prove (3.36) and (3.37). Note that by stationarity, (3.36) is 
a consequence of 

(3.38) lim \A T -A hom \=0, 

T— >oo 

for all d > 2, where £ • A T £, := ((£ + V0 T ) ■ A(£ + V^r)). Starting point for 

(3.38) is the definition of At and ^hom from which we deduce 

£ • (A T - A hom )C = <(£ + V0 T ) • A(£ + V</> T ) - (£ + V0) • A(£ + V0)> 

(3.39) = (£ • A(V</> T - V</>)> + (V<Ar • A(£ + V0 T )> 

_( V 0-A(£ + V0)). 

Let us treat each term separately. For the second term, we shall argue that 
(2.3) yields: for every stationary field £ : Z d — > R such that (£ 2 ) < oo, one has 

(3.40) T- 1 ^ T C) + (VC-A(e + V(/>T))=0, 

so that one may replace the second term of the right-hand side of (3.39) by 
— T~ (4>t) ■ For the first term, we shall use the following weak convergence of 
V(Pt{x) to V(j)(x) in probability: for every random variable % taking values 
in R rf with (|x| 2 ) < oo, one has for all x G Z d , 

(3.41) lim (x • (V<Mz) " V0(x))) = 0, 

T— >oc 

so that taking x = and x = A(0)^ shows that the first term in the right- 
hand side of (3.39) vanishes as Tfoo. For the last term, combining (3.41) 
and (3.40), we will prove 

(3.42) (V0-A(£ + V0)) = O. 

We directly draw the conclusion: the combination of (3.39), (3.41), (3.40) 
and (3.42) shows that 

lim sup |£ • (A T - A hom )£\ = lim sup T -1 ^), 

T->oo T->oo 

which implies (3.38) by Proposition 2.1. 

We give the arguments for (3.40), (3.41) and (3.42) for the reader's con- 
venience (we could also directly appeal to [13]). Multiplying the defining 
equation for 4>t by Q yields 

(3.43) T-^K) W - (V* • A(£ + Vcj )T ))(z)C(z) = 0. 
We then use the discrete Leibniz rule in the form 

V* • (CA(( + Vcp T ))(z) = (V* • A(C + Vcp T ))(z)((z) 

(3.44) 

d 

+ ^V*C(z)L4(£ + V0 T )(z-e J )] J . 
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Since V^t, £ and A are jointly stationary random fields, the expectation of 
the left-hand side of (3.44) vanishes, and 

((V* • A(t + V£r))(*)C(*)> = -fev*(( 2 )[i(( + VMz - ej))}\ 

(3.45) 

noting that V*jQ{z) = Vj({z — ej). We then take the expectation of (3.43) 
and use (3.45) to obtain (3.40). 

We recall the standard a priori estimate which one derives from (3.40): 

(T-Vt(x) 2 + |V^(:e)| 2 )<1. 

Since the left-hand side does not depend on x by stationarity, there exists 
g\% d — y~R. d such that up to extraction, V</>t(x) converges to g(x) weakly 
in probability for all x G Z d . By construction, g is a gradient field, and is 
jointly stationary with A. By the boundedness of (T -1 ^) 1 / 2 , one may pass 
to the limit in (3.40), and obtain for every stationary field £ 

(3.46) (VC-^(C + V0)) = O. 

As noticed by Kiinnemann in [13], this characterizes the gradient of the cor- 
rector, so that g = V</>. This proves (3.41) by definition of weak convergence 
in probability. 

We then use (3.46) for £ = 0^ and pass to the limit T f oo in (3.46) by 
the weak convergence (3.41). This proves (3.42). 

We finally turn to the proof of (3.37). By definition, (3.37) is proved if for 
all bounded random variables x> 

lim (x (£ + V0r(z)) • + V<£r(x))/i L (x) dx 

(3.47) 



X / (S + V<Kx))-A(x)(£ + V<l>(x))ii L (x)dx 

W.l.o.g. we may assume that x takes values in [0, 1]. By lower-semicontinuity 
of quadratic functionals in probability, and since x > 0, the weak convergence 
(3.41) of V^>t(x) to V0(x) in L 2 in probability for all x € Z rf yields 

liminf(x / (£ + V</> r (x)) • A(x)(£ + V</> T (x))/i L (x) dx 

fi L (x) (liminf + V^r(a:)) ■ A(s)(£ + V^r(ac))) J dx 



X I (f + V^(x)) • + V^(x)) ML (x) dx ). 
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Likewise 



liminf( (1 

T-s-oo \ 



x) [ {i + V(j> T {x))-A{x){i + V(j>T{x))ii L {x)dx 



) 




X) [ (€ + V<f>(x))-A(x)(£ + Vil>(x))ii L (x)dx 



since 1 — x ^ 0- Combined with the convergence of the expectation (3.36) 
and the trivial identity 1 = X + (1 — x), these two inequalities imply (3.47) 
for x taking values in [0, 1], and therefore (3.37) as desired. 

4. Proofs of the estimates on the Green functions. Before addressing 
the proofs proper, let us make a general comment. In what follows, we shall 
replace the classical Leibniz rule by its discrete counterpart. Although they 
are essentially the same, the expressions that appear are more intricate in 
the discrete case. In order to keep the proofs clear, we first present the 
arguments using the classical Leibniz rule (though it does not hold at the 
discrete level) and we later give a separate argument to show that the various 
results still hold with the true discrete version. 

4.1. Proof of Lemma 2.8. Without loss of generality, we may assume 
y = and suppress the y-dependance of Gt in our notation. We will first 
give the proof in the continuum case (i.e., using the classical Leibniz rule) 
and then sketch the modifications arising from the discreteness. 

We first argue that for any d, 



where for < M < oo, Gt,m denotes the following truncated version of Gt 

G T ,M = mm{G T ,M}>0. 
Indeed, we consider T~ 1 Gt — V* • AS/Gt = 5 in its weak form, that is, 



and select £ = Gt,m- Since Gt,mGt > G TM and provided that VGt,m • 
AVGt > VGt,m • AVGt,m, we obtain (4.1) by uniform ellipticity. Indeed, 
since A is diagonal, 



(4.1) 




(4.2) 




VG t ,m ■ AVG T (x) 



d 



= ^a(x + ei,x)(G T ,M(x + e { ) - G t ,m(x))(G t (x + e») - G T (x)) 
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> ^ a(x + ej, x)(G t ,m(x + e^) - G t ,m{x))' 



i=l 

2 



>a|VG T , M (x)| 2 . 
Step 1. Proof of (i) for d > 2. 

Following [7], Theorem 1.1, we argue that (4.1) implies a weak-L d// ( rf-2 ) 
estimate, that is, 

(4.3) £ d ({G T >M})<M- d ^ d - 2 \ 

For this purpose, we appeal to Sobolev's inequality in Z rf , that is, 

, , \ (d-2)/(2d) , r \ 1/2 

J zd G 2 T d Lf- 2) dxj <(jjVG T , M \ 2 dx) , 

which is a consequence of [22], Lemma 2.1 (when "re — t-oo"), or [3], Theo- 
rem 4.4 (when "r — > oo"). Via Chebyshev's inequality and (4.1), this yields 

MC d {{G T > M})^/^ < M 1 / 2 , 

which is (4.3). 

We now argue that the weak-L d ^ d ~ 2 ^ estimate (4.3) in Z d yields a strong 
L^-estimate on balls {|a;| < R} for 1 < q < More precisely, we have 

(4.4) f G q T dx<R d {R 2 ~ d ) q . 

J\x\<R 

Indeed, we have on the one hand 

r roc 

/ G q T dx = q C d ({G T > M'})M' 
Jg t >m Jm 

(4.5) + M q C d ({\G T \>M}) 

(4.3) 

< M q ~ d ^ d ~ 2 \ 

where we have used q < On the other hand, we have trivially 

(4.6) I G q T dx<R d M q . 

J{G T <M}n{\x\<R} 

With the choice of M = R 2 ~ d , the combination of (4.5) and (4.6) yields (4.4). 

In order to increase the exponent q in (4.4), one combines a Cacciopoli 
estimate 1 for monotone functions of Gt with a Poincare-Sobolev estimate to 
obtain a "reverse Holder" inequality (as in the proof of Harnack's inequality, 



,q - x dM' 



1 This is the only place where we use the Leibniz rule. 
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see [8], Chapter 4, Method II). We start with the Cacciopoli estimate, that 
is, 

(4.7) f \VG q J 2 \ 2 dx<R- 2 [ G q T dx 

J2R<\x\<m Jr<\x\<&r 

for all 1 < q < oo. For that purpose, we test (4.2) with £ = rj^GjT , where 
the spatial cut-off function r/ has the properties 

rj = l in {2R< \x\ <4R}, 

(4.8) 

?7 = outside {R<\x\ <8R}, iVr^iT 1 , 0<t?<1. 
This yields 

(4.9) T" 1 / 7] 2 G Q T dx+ [ V(r] 2 G q T ~ 1 )-AVG T dx = 0. 

Jz d Jz d 

Since by the uniform ellipticity of A, there exists a generic constant C < oo 
(which only depends on q, a, 0) such that 

V^GfT 1 ) -AVG T 

= {q- l)rj 2 Gf 2 VG T ■ AVG T + 2r]G q 7 l Vr] ■ AVG T 

Young 

> C^ri 2 G q f 2 \VG T \ 2 - CG q T \V V \ 2 
>C-W\VGf\ 2 -CG q T \V v \ 2 , 

[ T?\VG# 2 \ 2 dx< f G q T \Vr]\ 2 dx. 
Jz d Ji d 

In view of the properties (4.8) of rj, this yields (4.7) for d > 2. 
We now derive the "reverse Holder" inequality 

R~ d / G? /( ^ 2) dx) 

J2fl<|x|<4R / 

(4-10) " " yq 

<(R- d [ G q T dx) \ 

V JR<\x\<8R J 

For that purpose, we appeal to the Poincare-Sobolev estimate (see [22], 
Lemma 2.1, or [3], Theorem 4.4) on the annulus {2R < \x\ < 4R}: 

/ r N (d-2)/(2d) / , Nl/2 

( R~ d / \u\ 2d ^ d - 2 Ux ) < ( R 2 ~ d / \Vu\ 2 dx) 

V J2R<|z|<4fi / V J2R<\x\<4R / 

+ (R- d ( \u\ 2 dx^ 

J2R<\x\<4:R 



we obtain 
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4.3 



We apply the latter to u = G 9 / 2 : 

(d-2)/{qd) / r 

R- d I G «d/(d-2) d \ <(B?-d \VG q J 2 \ 2 dx 



2R<\x\<iR / V J2R<\x\<iR 



1/9 

(jrp djX J 

2_R<|x|<4_R 



+ {R- d I G q T dx 



The combination of this with (4.7) yields (4.10). 

We now may conclude in the case of d > 2: indeed, (4.10) allows us to 
iteratively increase the integrability q in multiplicative increments of 
in the estimate (4.4). Since any p < oo can be reached in finite multiplica- 
tive increments starting from a 1 < q < -r^, the side effect that the annuli 
get dyadically larger at every step does not matter qualitatively (in this 
sense, the above argument is much less subtle than the proof of the Harnack 
inequality). This proves (2.21). 

Step 2. Proof of (i) for d = 2. 

We now tackle the case of d = 2, which in fact amounts to the L 1 -BMO 
estimate 



(4.H) (V 2 / \u-u {{xl<R} \ q dx) ^ < [ 

V J\x\<R J J7? 



\f\dx 



\ 1/9 

u ~ v>t\x\<R\\ q di 

'\x\<R 

for 

(4.12) r 1 u-V*-iVu = /, 

where U{\ x \<r} denotes the average of u on the ball of radius R. We fix an 
exponent q < oo and a radius 1 R < oo and assume w.l.o.g. 

(4.13) H\x\<R} = °- 
As in (4.1), we have 

(4.14) [ \Vu M \ 2 dx<M [ \f\dx. 

J\x\<R Jl? 



As opposed to the case of d > 2, this is the only time we use the equation 
(4.12). 

Estimate (4.14) is used in connection with the Poincare-Sobolev inequal- 
ity with mean value zero, that is, 



r \ V s / r \ I/ 2 

R~ 2 I \um ~ (u M ){\ x \<R}\ S dx) < / \Vu M \ 2 dx 

J\x\<R / \J\x\<R 
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for any s < oo, which we use once for s = q, that is, 

(r~ 2 [ \um~ (u M ){\ x \<R}\ q dx) < ([ \Vu M \ 2 dx) 

\ J\x\<R J \J\x\<R J 

(4.15) 

(4.14) / [■ \ 1/2 



4.14) / r 

< i^M J \f\dx 



and once for arbitrary s (which we think of being larger than q) in the form 

l/s 

u M \ s dx ] 

*\x\<R 

1/2 



R- 2 [ \u M \ s dx] 

^ J\x\<R J 

(4.16) <([ \Vu M \ 2 dx\' +\(u M ){\x\<R}\ 

\J\x\<R / 
(4.14) / r \ V 2 / r \ 1/9 

< (m Jjftdxj + (^r 2 y 



'|s|<R 



/ /• \ s /( 2 <?) / /■ 

< m 1 ' 8 '^ ( \f\dx) + M 1_s / 9 ( RT 2 / 



We use (4.16) to estimate the peaks of u. More precisely, we claim that 
for s > 2q, 

1/9 

R~ 2 I \u\ q dx 

>{\x\<R}n{\u\>M} 

(4.17) 

s/(2q) ../_/• \ s /<? 2 

'M<-R 

The argument for (4.17) is similar to the case of d > 2: estimate (4.16) yields 
the weak estimate 

M{R~ 2 C 2 {{\x\ <R}n {\u\ > M})) l/s 

< (m [ \f\dx) ' +(r~ 2 [ \u\Ux) '\ 

\ Jl? J V J\x\<R J 

which we rewrite as 

R~ 2 £ 2 ({\x\ <R}n{\u\ >M}) 

(4.18) 

s/2 / r \ s /t 



\u\ q dx 
'\x\<R 



<M- s / 2 (f \f\dxY +M- s (r- 2 [ 

\Jl? J \ J\x 

On the other hand, we have 

p poo 

/ \u\ q dx = q C 2 ({\x\ < R}D{\u\ > M'})M' q ~ l dM' 

J {\x\<R}n{\u\>M} JM 

(4.19) 

+ M q c 2 ({\x\ < i}n{H >M}). 
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Since s > 2q, the combination of (4.18) and (4.19) yields 

RT 2 [ \u\ q dx 
J{\x\<R}n{\u\>M} 

/ r \ s/ 2 / r \ s /l 

< M qs/2 1 \f\dx) +M q - S (R- 2 \ U \ q dx 
\Jl? ) V J\x\<R 

which is (4.17). 

We now combine (4.15) and (4.17) as follows 

2 / \u\ q dx ) ^ 
'\x\<R J 

.13) / r \ 



(4.13) f f \ 

< [R~ 2 \u- (u M ){\ x \< R }\ q dx) 

\ J\x\<R J 

< \R~ 2 / \u M - (u M ){\ x \<R}\ q dx ) 

\ J\x\<R J 

+ (r- 2 [ \u\ q dx] ^ 

V J{\x\<R}n{\u\>M} / 

(4.15) and (4.17) / r \ 1/2 / r \ s/{2q) 

< Ml/2 U 2 i/i dx ) + Ml ~ s/{2q) U z2 i ; i dx 

M ls/q( R ~2 f \ u \1dx^ 



l\x\<R 

We claim that this estimate contains the desired estimate. Indeed, using the 
abbreviations 

U:=(r~ 2 [ {u^dx} 1 and F := / \f\dx, 
\ J\x\<R J JT? 

we rewrite the above as 

(4.20) U < M 1/2 F 1/2 + M 1_s/(2?) F s, ( 2q) + M 1 ~ s/q U s/q . 

Since s > q, choosing M ~ U sufficiently large, we may absorb the last term 
of (4.20) into the left-hand side. This yields 

U < U 1/2 F 1/2 + U 1 ~ s/ ^F s/ ^ 2q) . 

Using Young's inequality twice in the right-hand side since s > 2q, we obtain 
as desired U <F, which shows 

\G T -G Tlxl < R \ q dx) 1/q <l. 

\x\<R J 
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Step 3. Proof of (ii). 

We first derive a weak L 4 -estimate on \\x\ < R}: 

(4.21) ir 2 £ 2 ({G T > M} n{\x\< R}) < M~\ 
For that purpose, we combine (4.1), which for i?~ y/T turns into 

(4.22) R- 2 [ G 2 TM dx+ [ \X7G T , M \ 2 dx<M, 

with the Poincare-Sobolev estimate 

r \ 1/8 / f \ 1/2 

R- 2 \GTM-GTM SM <m\ 8 dx) < / \VG TM \ 2 dx 



M {\x\<R}\ 

\x\<R / \J\x\<R 

in form of 

, n 1/8 / , \l/2 

R- 2 / G 8 TjM dx) < / \VG T , M \ 2 dx) 

J\x\<R J \J\x\<R J 

R~ 2 I G 2 TM dx) ' . 

J\x\<R J 

This yields (4.21): 

(R- 2 M s C 2 ({G T >M}n{\x\ <R})) 1/8 < (r~ 2 [ G 8 TM dx) 1/8 ^M 1 / 2 . 

\ J\x\<R ' / 

We now argue that (4.21) yields (2.22). Indeed, combining 

R~ 2 G 2 T dx = qR~ 2 C 2 ({G T >M'}n{|x| <R})M'dM' 

J{G T >M}n{\x\<R} Jm 

+ R~ 2 M 2 C 2 {{Gt > M} n{\x\< R}) 



(431) 
< M 



with the trivial inequality 



R~ 2 G 2 TiM dx<M 2 

J\x\<R 



for M = 1 yields property (ii) of the lemma. 
Step 4. Proof of (iii). 
We establish for all q > 1 and R^> 1 

(4.23) (2R)~ d [ G q T dx<-^R- d [ G q T dx 

J\x\>2R. R JR<\x\<2R 
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Indeed, we test (4.2) with r/ 2 G^ _1 where the cut-off function 77 is chosen as 
follows 

77 = 1 in {\x\ > 2R}, 

rj = m{\x\<R}, |V?7|<-R _1 , < 77 < 1, 

yielding 

T~ l [ r] 2 G q T dx+ [ V(r ] 2 G q T - 1 )-AVGTdx = 0. 
Jz d Jz d 

Arguing as for (4.9), this yields 

T- 1 [ G q T dx+ [ \VG q rl 2 \ 2 dx<ir 2 [ G q T dx, 

J\x\>2R J\x\>2R JR<\x\<2R 

so in particular (4.23). 

We now turn to (2.23). We introduce the abbreviations 

R k := 2 k VT, 



A k :=R~ d G q T dx, 
JR k <\x\<R k+1 



so that (4.23) turns into 



T 



Afc+i < C-j^ A k = C4 k A k , 
K k 

where C denotes a constant depending only on a, (3, and d. This yields by 
iteration 

fc-i 

Afc < AoC fc Yl 4 ~* = A C k 4:-( k - 1)k/2 = AqC^-^- 1 ^. 
Thus, for all r > 0, 



In I (-7=) J < A;rln2+(fe + l)lnC-^ln2 
< 1 

for k large enough. Hence, 

/ G q T dx<A Ri(^) r . 

To conclude the proof of (iii) , it remains to argue that 
(4-24) A <{j = | \^ d)q _ 
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For d > 2, this a consequence of (2.21), whereas for d = 2 we combine (2.20) 
with (2.22) as follows: 



A < T- 1 [ G T (x) q dx 

J\x\<2VT 

£r '(((^ 1 "- J w*) 



(2.20) 

< Vt 2 



t^, (Lte ) ' J ') , ") , 



(2.22) 

< Vt 2 



< 1. 



5iep 5. Modifications due to the discreteness. 

The only place where we have used the Leibniz rule is the proof of the 
Cacciopoli inequality (4.7). At the discrete level, we have for £ € {1, ... , d} 

V^Gf 1 )^) 

= V 2 (x + ei)G g T -\x + e;) - if (x)G q ~ l (x) 

(4.25) 

= riix + e^+yjx) ^^ + ^ _ Qq - 1{x)) 

+ 7?2(x + e ; ) ~ 7?2( " ) (^- 1 (^ + *) + Gf\x)). 

Taking advantage of the diagonal structure of A (although this is not essen- 
tial), we obtain 

d 

V( V 2 G q f l )-AVG T (x) = ^VitfGZr^Wafax + edViGrix) 

i=i 

(4.25) A 7] 2 (x + ei) + rj 2 (x) 

= 2^ a ( :E ' ;E + e i) 2 

1=1 

x (Gfr 1 (x + e l )-G«r 1 (x))V 4 G T (x) 



>0 

r/ 2 (a; + e^) — rj 2 (x) 

i=i 



+ y^a(x,x + e$)- 
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x (G g T -\x + ei) + G«T 1 (x))V l G T (x). 

Since the underbraced term is nonnegative, the lower and upper bounds on 
a yield 

V(? ? 2 G^ 1 ) • AVGt(x) 

> a £ T?2(X + e ^ + ??2(x) (Gfr 1 (^ + *) " GT 1 (^))V l G T (x) 

i=l 

</ 



- ^ IV^) ^ + ^ + V{X) (G«-\x + e,) + G$r\xWiGr(x)\ 

i=l 

> a 2^ 2 ^ ^ (x + ejj-Gf (x))ViG r T (a;) 

i=l 

- /3C^(G T (x + + G T {xY)\V iV {x)\ 2 

i=i 

-^g( ^ + e - ) + ^ ) ) 2 (V,G r (^)) 2 

i=l v ^ , 

<(r,2(x+ei)+ V 2(x))/2 

x (Gfr 2 (x + ei ) + G^ 2 (2;)). 
Using the inequality (proved at the end of the step) 

2(69-1 _ c q-i^ b _ c ) > (5 _ c f{b q ~ 2 + c^ 2 ) 

(4.26) 

for b,c>0,q>2, 

we may absorb the last term of the right-hand side of the latter inequality 
into the first term for C large enough, so that it turns into 

VtfGf 1 )- AVG T {x) 

>(q-2/3C- 1 )X: 7?2(x + e - ) + 7?2(x) 

(4.27) 



x (GfT i (x + e l )-Gfr i (x))V l G T ( 

d 

- f3C^2(G T (x + eO 9 + Gt^^IV^x)! 2 . 

8=1 

Using now the following inequality 

(4.28) (6 g - 1 -c 9 - 1 )(6-c)>(& ,?/2 -c 9/2 ) 2 for b,c>0,q>l, 



x 
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(4.27) finally turns into 

V(rj 2 Gf l ) -AVGt(x) 



> J2 v2{x + e ^ +7]2{x) (Gf(x + *) - Gf( X )Y 



2 
=i 

d 



- CY,(Gt(x + e t y + G T {xY)\V iV {x)\ 2 . 

i=l 

Combining this with (4.9) yields 

r) 2 (x)\VG^ 2 (x)\ 2 dx 

[ (G T (x + ei y + G T (x) q )\V iV (x)\ 2 dx, 
Jz d 



(4.29) 

< 



which implies as desired 

f \VG q J\x)\ 2 dx<R- 2 f G T (x) q dx, 

J2R<\x\<4:R JR<\x\<8R 

provided that rj satisfies in addition 

rj(x) = for xi {y.R+ 1 < \y\ < 8R- 1}, 

which is no restriction since R S> 1 • 

We quickly sketch the proofs of (4.26) and (4.28) to conclude. Inequality 
(4.26) follows by symmetry from 

( 6 <z-i_ c <?-i)(5_ c )_( 6 _ c )2 cg -2 

= (b-c)(b q - 1 -bc q ~ 2 ) 
= b{b-c)(b q ~ 2 -c q ~ 2 ) 
= b\b-c\\b q - 2 -c q ~ 2 \ >0. 

To prove (4.28) we first note that by homogeneity and nonnegativity of b 
and c, it is enough to consider c = 1 and b > 0. We introduce the function 
h = R+ ->• R+ defined by 



h(b) = < 



U(«-i)' 
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Since h > 0, the claim is proved if h is bounded on R + . As /i(0) = 1 and 
lini;,_s l00 h(b) = 1, it is enough to prove that h is continuous on R + . A Taylor 
expansion around 6 = 1 yields 

(& 9/2 -l) 2 = ^-l) 2 + o((6-l) 2 ), 
(69-1 _ i)( 6 - 1) = ( q - i)( 6 - 1)2 + (( 6 - 1)2). 

Hence, lim^i h(b) = h(l), h is continuous and therefore bounded on R + , as 
desired. 

4.2. Proof of Lemma 2.9. The proof relies on three ingredients: a Mey- 
ers' estimate based on the L q theory for the constant-coefficients Helmholtz 
projection, a Cacciopoli estimate and the estimates of Lemma 2.8. 

We begin with Meyers' estimates. Let u : Z rf — > R, / : Z d — > R, and g : l* d — > 
R d have support in {\x\ < R}, and let u satisfy the equation 

(4.30) -V* ■ A{x)Vu(x) = V* ■ g{x) + f(x) inZ d . 

We claim that there exists p > 2 depending only on a,/3, and d such that 
for all R S> 1 , the following L p -estimate holds 

\ 1/p / /■ \ !/p 

|V«(x)| p da; < / |g(x)| p dx 



(4.31) + (7 |/(x)| 2 dx) 

As in the original paper [16] by Meyers, the proof of (4.31) relies on a 
perturbation argument and on the L q regularity theory for the Helmholtz 
projection. 

Step 1. L q regularity for the Helmholtz projection. 

Let V.:L 2 (Z d ,R d ) -)■ L 2 (Z d , R d ) denote the Helmholtz projection, that 
is, the orthogonal projection onto gradient fields for the inner product of 
L 2 (Z d ,R d ). By definition, U is continuous on L 2 (Z d ,R d ) and satisfies 

( 4 - 32 ) ll^9llz,2(Z d ,R d ) - lbllL2( Z d )R d). 

Let us show that T~L can be extended to a continuous operator from L q (Z d , R d ) 
to L q (Z d , R d ) for all 1 < q < oo. The proof is standard, appealing to Calderon- 
Zygmund singular integral theory and to Marcinkiewicz interpolation theo- 
rem (such theorems apply to the discrete case under investigation since the 
associated measure has the so-called "doubling" property). Since H com- 
mutes with translations, it is a convolution operator: there exists a matrix- 
valued kernel K such that 

(4.33) Hg{x)= I K(x - y)g{y) dy. 
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From an elementary Fourier series analysis (see [15] for related arguments), 
we infer that the symbol of K coincides with the symbol of the second 
derivative of the Green's function of the Laplace equation studied in [15]. In 
particular, from the analysis of [15], we learn that 

( 4 - 34 ) |Vif(x)|< 1 

1 + |x| a+i 

We are therefore in position to apply Calderon-Zygmund's theory (see [20], 
Theorem 2, page 17), which shows that T~L is of weak type (1,1) (see the proof 
of [20], Theorem 3, page 19). Appealing to Marcinkiewicz' interpolation 
theorem (see [1], Theorem 1.3.1, page 9) then shows that "H can be extended 
to a continuous operator from L q (Z d ,R d ) to L q (Z d ,R d ) for all 1 < q < 2. A 
standard duality argument (see, [19], 2.5(c), page 33, e.g.) implies that H can 
also be extended to a continuous operator from L q (Z d ,R d ) to L q (Z d ,R d ) for 
all 2 < q < oo. Let r > 2 be fixed, and for all q > 1 let denote by C q the norm 
of U in C(L q (Z d , R d ),L q (Z d , R d )). Then Riesz-Thorin interpolation theorem 
(see [1], Theorem 1.1.1, page 2) shows that for all 6 G (0,1), C 2 e+ r (i-e) < 
C e 2 C}r 9 , so that 

(4.35) limsupC 3 <l 

<j-s>2 

since C 2 < 1 by (4.32). 

We now turn to the proof of (4.31) proper and proceed with the pertur- 
bation argument. 

Step 2. Proof of (4.31) for / = 0. 

We first assume that / = 0, and rewrite the left-hand side of (4.30) as a 



perturbation of the operator — ^j^A: 



2 ' 

a + /3 An = V*.( 5+ (^-^Id)v.) 
or equivalently in the form 

(4.36) _A. = V.(^(,+ (A-H±^ld)v.)). 

In order to apply the L q theory for the Helmholtz projection, we need to 
show that 

2 ( t. a+£ V 



(4.37) Vu^—[3 +[ A-^Ii)Vu 

Since Vm is obviously a gradient, it remains to show that for all £:Z d 
such that VC G L 2 (Z d ,R d ) one has 

/ Vu(x) ■ V((x) dx 

Jw,d 
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(4.38) 



a + /3 



9+[A 



ot + 



Vu (x) • VC(x) dx. 



To this aim, we multiply (4.36) by £ and integrate by parts using that u, 
Vu and g have compact supports. This yields (4.38) and proves therefore 
(4.37). The continuity of U from L q (Z d ,R d ) to L q (Z d ,R d ) proved in Step 1 
then implies that 



(4.39) 



\x\<R 



\Vu(x)\ q dx 



1/9 



a + 13 



\x\<R 



^^t^ldjVu(x) 



9 \ V« 



Using the triangle inequality, (4.39) turns into 
\Vu(x)\ q dx 

\x\<R ) 



(4.40) 



a + /3 



\ 1/9 

b(x)|«dx 
|x|<i? / 



" + V|x|< 



A(x) 



Since a 6 -4, Q/ g, | (j4(x) — ^^Id)Vu(x)| < ^2^|V«(x)| and we may absorb 
the term 



a + /3 



§-a\ 



Id ]Vu(x) 



9 \ V? 

dx 



a + /3 



|ac|<J2 

/9-a 



A(x) 



a + /3 



Id 



]Vu(x) 



9 \ 1/9 

fix J 



|Vu(x)| 9 (ix 

« + /3V4|<R 
into the left-hand side of (4.40) provided that 



1/9 



(4.41) 



c p ~ a CI 

<i 



The interpolation property (4.35) ensures there exists p > 2 such that (4.41) 
holds for all p > g > 2. For such & q, we then have 



(4.42) 



1/9 

\Vu(x)\ q dx\ < 



\x\<R 



\x\<R 



\g(x)\ q dx 



1/9 
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as desired. 

Step 3. Proof of (4.31) for general /. 

Note that since u and g have compact supports, equation (4.30) implies 
that f zd f(x) dx = 0. We first show that there exists Vu> G L 2 (Z d ,R d ) such 
that for all ( : % d -> R with VC G L 2 (Z d , M d ), one has 



(4.43) / Vw(x) ■ V((x)dx= / f(x)((x)dx, 

Jz d Jz d 

so that (4.37) turns into 



V. = M^^+^--^ia)V«) + V,. 
Provided 

(4.44) ^ |Vt«(s)|«daA < R 1 -** 1 /*- 1 /*) (J fixfdx^' , 

for all 2 < g < g for some q> 2, we then conclude as in the case / = (with 
potentially a smaller p). To prove the existence of such a Vu;, we proceed 
by minimization and consider the problem 



inN / \V((x)\ 2 dx- / f(x)((x)dx; 
[Jz d Jz d 

(4.45) 

C:Z d ^R,V(€ L 2 (Z d ,R d ) 



The same argument as in the proof of Riesz' theorem yields the existence 
of a minimizer once one shows that the functional is coercive. Let R be 
large enough so that / has support in {|x| < R}, and denote by C{|x|<i?} t ne 
average of C on {|x| < R}. Since / has zero average, one may subtract the 
average of C and obtain by Cauchy-Schwarz and Poincare's inequalities 

f(x)((x) dx 



f(x)(((x) - C{\ x \<R})dx 

\x\<R 

<r( [ f(x) 2 dx) ' ( [ \V((x)\ 2 dx X 



\x\<R 



<-2R 2 [ f{x) 2 dx + \ [ |VC( 

JZ d z JZ d 

This shows that for all test functions £ 

\V((x)\ 2 dx- [ f(x)C(x)dx 



x)\ 2 dx. 



(4.46) 

> 



-2R 2 [ f( x ) 2 dx + l [ \V((x)\ 2 dx, 
JZ d 1 JZ d 
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as desired. This proves the existence of a minimizer w : Z d — > R such that 
Vu; G L 2 (R d ,Z rf ). In addition, it satisfies the estimate 

(4.47) f \Vw(x)\ 2 dx<4R 2 [ f(x) 2 dx. 

Jz d Jz d 

Since w is a minimizer of (4.45), the first variation of the energy at w van- 
ishes, and w satisfies (4.43). 

It remains to estimate the L q norm of Vw for some q > 2. To this aim, 
we argue that 



\VVw(x)\ z dx = I (Aw(x)) z dx 

(4.48) 

f(x) 2 dx. 



As in the continuum case, the first identity in (4.48) follows directly from 
two integrations by parts for w with compact support. For general w, the 
boundary term involves products of first and second derivatives of w on 
spheres of large radius R. In our discrete setting, these boundary terms can 
be estimated by the integral of |Vw| 2 outside the ball of radius R, which is 
finite since |Vu;| 2 is integrable by construction. Hence, the boundary terms 
can be made to vanish in the limit R — > oo. The second identity in (4.48) 
follows from the fact that w solves the equation 

-Aw{x) = f(x) inZ d , 

which is a consequence of (4.43). We are in position to conclude. For d > 2, 
we appeal to Poincare-Sobolev inequality on Vw to turn (4.48) into 

(4.49) / \Vw(x)\ 2d ^ d - 2 Ux< [ f(x) 2 dx. 



Combined with (4.47), (4.49) implies (4.44) for all 2 < q < ^ by Holder's 
inequality. For d = 2, we appeal to Poincare-Sobolev inequality on (V^u;) 2 
for i G {1, . .. ,d} to turn (4.47) and (4.48) into 

/ (V l w(x)) 4 dx<( [ iViiVM^fldx] 

(4.50) < / \VVw{x)\ 2 dx / \Vw(x)\ 2 dx 

Jz 2 Jz 2 



<R 2 ^fj{x) 2 dx^ 



2 



Combined with (4.47), (4.50) implies (4.44) for all 2 < q < 4 by Holder's 
inequality. 
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Step 4. Cacciopoli estimate. 

We need the following finer version of (4.7): For all k£R, 

\VG T (x)\ 2 dx 

2R<\x\<16R 

(4.51) <R~ 2 [ (G T {x)-K) 2 dx 

J R<\x\<32R 

+ T~ 1 \k\ / \G T (x)-n\dx. 

JR<\x\<32R 

This variant of Cacciopoli's estimate can be proved along the lines of (4.7), 
multiplying the equation by i] 2 (Gt(x) — k) instead of t] 2 Gt(x). The zero 
order term then brings the new term in the right-hand side of (4.51). By 
Young and Cauchy-Schwarz' inequalities, the second term of the right-hand 
side is controlled by 



T 1 \k\ / \Gt{x) — k\ dx 

JR<\x\<32R 

= T- 1 \k\R [ RT l \G T {x)- n\dx 

J R<\x\<32R 

<T- 2 \k\ 2 R 2 + R' 2 [ \G T {x)- K \ 2 dx. 

Jr<\x\<32R. 



lR<\x\<32R 

Hence, it only remains to estimate the first term of the right-hand side 
of (4.51). To this aim, we appeal to Holder's inequality with exponents 
(p/2,p/(p-2)) forp>2: 

f \G T (x)-K\ 2 dx 

JR<\x\<32R 

<( f \G T (x) - K\ p dx) ' P (R d )tP-V/p 

\jR<\x\<32R J 

G T {x)-K\"da 

lR<\x\<32R 

Using these last two estimates and the elementary inequality (a 2 + 6 2 ) 1 / 2 < 



R d ~ 



~ 2d /p( [ \G T (x) - k\ p dx) 

\Jr<\x\<32R J 

t 

a + b, (4.51) turns into 



1/2 

\X7G T (x)\ 2 dx 

2R<\x\<16R 



(4.52) <R- 1 R d ' 2 - d / p ( f \G T {x) - k\p dx) ' 

\jR<\x\<32R J 
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+ T~ 1 R d / 2+1 \K\, 

that we will use with k = Gr{R<|a;|<32_R}- 

Step 5. In this step, we use Steps 1 and 4 to argue that 

i/p 

\VG T (x)\ p dx' 

4R<\x\<8R 



(4.53) <(J?" 1 +l£r- 1 )f / \G T (x)-G T{R < ]x] < 32R} \ p dx) 

\J2R<\x\<16R J 



\ VP 

Gt(x) - Gt{R<\ x \<32R}\ 1 ' ds 

l2R<\x\<16R 
x\<32R}- 

We apply Meyers' estimate (4.31) to the function u = ti(Gt — Gt {r<\ x \<32R}) > 
where the cut-off function r\ : Z, d — > [0, 1] is such that 

r}(x) = l for AR< \x\ < 8R, 

(4.54) 

^)=° io *{t\iwR-i, w*^. 

For all i G {1, . . . , d}, the discrete Leibniz rule yields 

Vi«(x) = r)(x)ViG T (x) + (G T (x + efi - G T { R <\ x \< mR })ViT]{x). 
Based on this, a direct calculation shows 
-V* • AVu(x) 

d 

= - r](x)V* ■ AVG T (x) - y^V*r](x)a(x-e i ,x)V*GT(x) 

(4.54) and (2.11) 1=1 
= ri(x)T L G T (x) 

d 

~ V *i (( Gt ( x + e «) ~ G T {R<\ x \< 3 2R})a(x, X + ei)Vi7](x)) 
1=1 

= V* • ^-^(G^(s + e i )-G r T{fl<|x|<32fl})a(«>a ; + ei)Vi»7(s)e i ^ 

d 

- Yl V*iV(x)a(x - e u x)V*G T {x) - r ? (x)T^ 1 G r (x). 

i=l 

This identity has the form of (4.30) provided we define the functions / and 
9 by 

d 

f(x) = -Y,^Mx)a(x - e u x)V*G T (x) - ^T^G^x), 

i=l 
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d 

9{x) = -^2(G T (x + e { ) - G T { R <\ x \<32R})a(x,x + e i )V i r](x)e i . 
i=i 

Since u, f and g have support in {\x\ < 16R}, we may apply estimate (4.31) 
which yields 

\ i/j> 

\Vu{x)\ p dx 

'|x|<16R / 

\ VP 

\x\<32R}\ P dx 



£ fe / WiV{x)\ p \G T {x + e,)-G T{R ^ 
\ i=i ~' zd 



+ R l-d(X/2-X/p) 

r \ 1/2 

/ (\V* V (x)\ 2 \V*G T (x)\ 2 + T- 2 r](x) 2 G T {x) 2 )dx) . 
Jz d ) 

Using the property (4.54) of rj, and the triangle inequality, we are left with 

\ i/p 
\Vu(x)\ p dx) 

'|a;|<16fl / 

X/p 



a 

< W / \Gt{x) - G T{R < lxl < 32R} \ p dx) 

\J2R<\x\<X6R / 

(4.55) 

/ r x 1 / 2 

+ R d/ P -d/2l / \VG T (x)\ 2 dx) 

\J2R<\x\<X6R J 

+ R d/ P -d/2+if f T- 2 G T {x) 2 dxX' 

\J2R<\x\<Xm J 



l2R<\x\<X6R. 

Let us rearrange the terms. For the third term, the triangle inequality and 
Holder's inequality with exponents (p/2,p/(p — 2)) show that 

(( G T {x) 2 dx\' 

\J 2R<\x\<X&R, J 
R^ 2 Gt{R<\ x \<32R} 

+ R d/2-d/pf f \G T (x)-G T{R<lxl<3 2R } \ P dx) 1/ \ 
\J2R<\x\<X6R J 

whereas for the second term we appeal to the Cacciopoli estimate (4.52) 
with k = Gt{r<\ x \<32R}- Hence, (4.55) finally turns into 

X/p 

\X7u(x)\ p dx ' 

\x\<X6R 
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(4.56) <(i?- 1 + i?T- 1 )f / \G T (x)-G T{R < lxl < 32R} \ p dx) /P 

\J2R<\x\<16R J 
+ T~ 1 R d / p+1 G T {R<\ x \<32R}- 

We are in position to conclude the proof of this step. For all i € {1, . • • ,d}, 
the discrete Leibniz rule yields Vjtt(x) = rj(x)S7 \Gt{x) + Gt(x + ei)Virj(x). 
Hence, (4.54) implies that Vu(x) = VG T (x) for 4R < \x\ < 8R, so that (4.56) 
yields (4.53). 

Step 6. Proof of (2.24). 

We claim that (2.24) follows from (4.53) and the estimates of Lemma 2.8. 

We distinguish two regimes: R < y/T and R > y/T. We begin with R < 
y/T. For the first term of the right-hand side of (4.53), we appeal to the 
BMO estimate (2.20) of Lemma 2.8 for d = 2 and to the decay estimate 
(2.21) with "q =p" for d > 2, so that 

{R- l + RT~ l )(f \G T (x)-G T{R < lxl < 32R} \ p dx) ^ 

\J2R<\x\<16R J 

(4.57) " " 

< R-\R d R( 2 - d )P) 1/p = R d /p- d +\ 

For the second term, we estimate the average using (2.22) for d = 2 

_ (2.22) 
G T{R<\x\<32R}^R~ 2 ^T 2 Gt { \ x \< 32 ^t} ^ R~ 2 T, 

and using (2.21) with u q = 1" for d > 2 

— < ' 2 '' 21 ^ 2d d 

Gt{R<\ x \<32R} ^ R <R 

since R < y/T. Hence, in both cases, 

(4.58) T^R^Gt^^r} < R d /P- d+1 . 

From (4.57) and (4.58), we then conclude that (2.24) holds for R < y/T. 

We now deal with the case R > y/T. For the first term of the right-hand 
side of (4.53), we use the decay estimate (2.23) with exponents "q =p,r = 
k + 2p" which yields 

/ \G T (x) - G T {R<\ X \<32R}\ P dx < / G T (x) p dx 

J2R<\x\<16R JR<\x\<32R 

(2 23) 

< R d R^ p (VTR- 1 ) k+2p , 

and therefore 

{R~ l + RT~ l )(f \G T (x)-G T{ R<\x\<32R}\ P dx) ? 

\J2R<\x\<16R J 

(4.59) 

<R d / p - d+1 (y/TR- 1 ) k/p . 
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For the second term, we proceed the same way, and appeal to (2.23) with 
exponents u q = 1, r = k/p + 2," which yields 

G T{R < lxl <32R} < R^iVTR- 1 )^ 2 = TR^VTR^, 
and therefore 

(4.60) T- l R d /P +1 G T{R ^ 32R} < R^+WTR- 1 )^. 

From (4.59) and (4.60), we then deduce that (2.24) holds for R > VT as 
well. 

4.3. Proof of Corollaries 2.2 and 2.3. These results are easy consequences 
of Lemmas 2.8 and 2.9. We include their proofs for convenience. 

4.3.1. Proof of Corollary 2.2. W.l.o.g. we assume y = and skip the 
dependence on y in the notation. We distinguish two regimes: \x\ < VT and 

\x\ > Vt. 

In the first case, we use (2.22) and the intermediate results (4.4) in the 
proof of Lemma 2.8, which yield 



for d 



2 / G 2 T (x)dx<T, 
J\x\<Vt 

ford>2 / G q T (x)dx< VT d (VT' 

J\x\<VT 



i\x\<Vt 

and imply for q = j=± G (1, ^) by the L 2 - L°° estimate 

(4.61) G T (x) < VT for \x\ < VT. 

For \x\ > VT, we use the decay estimate (2.23) of Lemma 2.8 with u q 
d,r = d(d + l) + l" 

G d T {x) dx < R d {VTR- 1 ) d{d+1)+1 = VT d {VTR- 1 f +1 , 



L 



<R<\x\<2R 

so that we may deduce 

(4.62) G T {x) < (VTR-if^VT for R < \x\ < 2R. 

We then define h T G L l (R d ) by 

VT2- k{ - d+1 ' d \ VT2 k <\x\< VT2 k+1 , ken, 
VT, \x\ < VT, 



h T (x) 



so that Gt(x) < hx(x) for all x G Z d . This concludes the proof since the 
factors in (4.61) and (4.62) only depend on a,f3 and d. 
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4.3.2. Proof of Corollary 2.3. We divide the proof in three steps. We first 
prove that the Green function Gt{x,u) is symmetric so that V x GT{x,y) = 
V x GT(y,x). In the second step, we show the uniform bound for \x — y\ > R 
sufficiently large, and in the third step for \x — y\ < R. 

Step 1. Symmetry of Gt- 

Let y,y S Z d . Testing the defining equation (2.11) with x h-> GT{x,y) 
yields 

f T- 1 G T (x,y)G T (x,y)dx+ [ VG T (x,y) ■ A(x)VG T (x,y)dx = G T (y,y). 

Since A is symmetric, the left-hand side of this identity is symmetric in y and 
y. Hence, the right-hand side is also symmetric, that is, Gx(y, y) = Gt(v, y)- 

Let R ~ 1 be sufficiently large so that Lemma 2.9 applies. 

Step 2. Estimate for \x — y\ > R. 

For q = 2, formula (2.24) yields for all k € N 

/ \V x G T (x,y)\ 2 dx<(2 k R) d ((2 k R) 1 - d ) 2 = (2 k R) 2 - d < 1. 

J2 k R<\x-y\<2 k + 1 B. 

Hence, by the discrete L 2 — L°° estimate, this shows 

(4.63) \V x G T (x,y)\<l for \x - y\ > R. 
Step 3. Estimate for \x — y\ < R. 

We now use an a priori estimate. Let i € {1, . . . , d} be fixed. We set u(x) := 
Gt(x, y + ej) — Gt(x, y) = ^yiGr{x, y). This function solves the equation 

(4.64) T~ 1 u - V* • AVu = f mZ d , 

where f(x) = 5{y + — x) — 5(y — x). Since / satisfies f zd f{x) = 0, one has 
by integration by parts, elhpticity of A and Poincare's inequality 



/ T 1 u(x) 2 dx + a / \\7u(x)\ 2 dx < / f{x)u(x)dx 
J% d Jz d Ji d 

f{x){u(x)-u { \ x _ y \< R] )dx 
<r([ \Vu(x)\ 2 dx) ' . 

\J{\x-y\<R} J 



L 



'{\x-y\<R} 

Hence, 

\Vu(x)\ 2 dx< R 2 ~l. 



This shows that sup|Vu| < 1. Therefore, for all x such that \x — y\ < R, we 
have using Step 2 and the fact that R is of order 1 

\u(x)\ < i?sup| Vn| + sup |n(z)| < 1. 

| 2 -2/|>R 



G2 
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Recalling that u(x) = V Vi GT(x,y), we conclude by Step 1 that this implies 
\V yi G T (y,x)\ < 1, as desired. 

5. Proofs of the other auxiliary lemmas. 

5.1. Proof of Lemma 2.3. W.l.o.g. we may assume 



(5.1) 



E 

i=l 



sup 



dX 



dcii 



< oo. 



Let X n denote the expected value of X conditioned on a\, . . . ,a n , that is, 

X n (ai,. . . ,a n ) = (X\ai,...,a n ). 
We will establish the following two inequalities for n < n G N: 



(5.2) 



(5.3) 



(X 2 n ) - (X n ) 2 < J2 



i=l 



sup 



dX 



da; 



((Xn - X n ) 2 ) < ( sup 

i=n+l > a ' 



dX 



dai 



var[ai], 



var[ai]. 



Before proving (5.2) and (5.3), we draw the conclusion. There is a slight 
technical difficulty due to the fact that there are infinitely many random 
variables. 

From (5.3) and (5.1), we learn that {^nlntoo is a Cauchy sequence in 1? 
w.r.t. probability. Hence, there exists a square integrable function X of a 
such that 



(5.4) 



lim((X-X n ) 2 >=0. 

ntoo 



By construction of X n , (5.4) implies 

(X\a±, . . . , a n ) = (X\ai, . . . , a n ) for a. e. (a\, . . . , a n ) and all n £ N. 

This means that the random variables X and X agree on all measur- 
able finite rectangular cylindrical sets, that is, measurable sets of the form 
A\ x • • • x A n x R x • • • , where n is finite. Since these sets are stable under 
intersection and generate the entire cr-algebra of measurable sets, the ran- 
dom variables X and X are uniquely determined by their value on these 
sets [10], Satz 14.12. Hence, the two random variables coincide, yielding 

(5.5) X = X almost surely. 

From (5.2), (5.4) and (5.5), we obtain in the limit n|oo as desired 



var[X] = (X 2 ) - (X) 2 < 



oo 

E 

i=l 



sup 



dX 



da,. 



var[ai]. 



VARIANCE ESTIMATE FOR EFFECTIVE DIFFUSION 63 
We now turn to (5.2) and (5.3). Notice that we have the decomposition 

n 



i=l 



where we have set Xq : = {X) so that {X n ) 2 = {Xq). Hence, (5.2) reduces to 



(5.6) 

Likewise, 



(Xt)-(XU)< sup 



dX 



var[ai]. 



{{x n - X n f) = (xl) - {xl) = «x?> - {xl x )), 



i=n+l 



so that also (5.3) reduces to (5.6). 

We finally turn to (5.6). We note that by our assumption that {oj}j S N are 
i.i.d., we have 

(X?(ai,...,ai)) = (^j Xf(a 1 ,...,a i - 1 ,a^P(da' i )\ 

Xi-\{ax,...,ai-i) = j Xi(a!,. . . ,ai-i,a")P(da"), 
where (3 denotes the distribution of a±. Hence, we obtain 

(X?) - (xl,) 

= (/ X i( ai > ' " ,0 ' i - 1,a ' i )P( da ' i ) ~ (/ X i( a U---> a i-U a i)P( da i)) \ 

^(Xi(ai, . . .,ai-i,a'i) - X^m, . . . ,ai_i,a")) 2 ^(da-)/3(da-')\ 



< 



sup 



sup 

a!" 



sup 

< 

dXi 



dXi 



da; 



(ai,...,a,_i,aj ) 



dcii 
dXi 



den 



(oi, . . . , ai-i,a'l') 



(oi, . . . , ai-i,a'l') 



'^[-allf^da'mda'l) 



(a-) 2 /3(da-) - ( / a"/3(da" 



var[oJ. 



We conclude by noting that by the definition of Xi and Jensen's inequality 



dXi 



da; 



(ai,...,Oi) 



j dX 
\ dai 



Ol , . . . , Oj 



< 



dX 



dai 



ax,...,a t I , 
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so that 
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sup 



(at, . . . } Oi_i,<4) 



< ( ( sup 



( sup 

\ a! 



dX 

da; 



Oj-i, a'n a>i+i, ■ ■ •) 



(ai, . . . ,0,-1,^,0^+1, . . .) 



ai , . . . , a, 



5.2. Proof of Lemma 2.5. Let us divide the proof in four steps. 

Step 1. Proof of (2.15). 

We recall the definition of the operator 

(Lu)(x) = a(x,x')(u(x) — u(x')). 

x',\x'— x\=l 

For convenience, we set e = [z, z'], z' = z + ei. We recall that Gt(-,v), y £ %> d , 
is defined via 

(5.7) (T- 1 + L)G T (;y)(x)=5(x-y), x £ Z, d . 

Hence, we obtain by differentiating (5.7) 

( (T_1 + L) dajej Gt( '' y) ) (X) + y) " Gt(/ ' y))< ^ " Z) 

+ (G T (z',y)-G T (z,y))8(x-z') = 0, 
which, in view of (5.7), can be rewritten as 
d 



(5.8) 



(T- 1 + L) ( -_G T (-,y) + (G T (z,y) - G T (z',y))G T (;z) 



da(e) 



+ (G T (z / ,y)-G T (z,y))G T (;z') ) j = 0. 



From this, we would like to conclude 



(5.9) 



<9a(e) 



G T (;y) + (G T (z, y) - G T (z', y))G T (-, z) 



+ (G T (z',y)-G T (z,y))G T (;z') = 0, 



which is nothing but (2.15). 

In order to draw this conclusion, we will appeal to the following uniqueness 
result in L 2 (Z d ): any u £ L 2 (Z d ) which satisfies ((T _1 + L)u){x) = for all 
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x G TL d vanishes identically. However, we cannot apply this directly to u given 
by the left-hand side of (5.9), since we do not know a priori that gx^ T^)!/) 

is in L 2 (Z d ). 

For that purpose, we replace the derivative -q§^ by the difference quo- 
tient. We thus fix a step size h 7^ and introduce the abbreviations 

G T {x,y) :=G T (x,y;a) and G' T (x,y) := G T {x,y,a), 

where the coefficients a' are defined by modifying a only at edge e by the 
increment h, that is, 

a'(e) = a(e) + h and a(e') = a(e) for all e 7^ e. 

We further denote by Lj- :=T + L a and L' T := T + L a i the operators with 
coefficients a and a', respectively. We mimic the derivation of (5.8) on the 
discrete level: from (5.7), we obtain 

= ^L T G T (;y)-L' T G' T (;y)) 

= L T ±-(G T (;y) - G' T {-,y)) + \(L T - L' T )G' T (;y) 

= L T ±-(G T (-,y) - G' T (-,y)) + (G' T (z, y) - G> T (z> ,y))5(- - z) 
+ (G' T (z',y)-G' T (z,y))5(--z') 

= L T (±(G T (;y) - G' T (;y)) + (G' T (z,y) - G' T (z' ,y))G T (;z) 

+ {G' T (z',y)-G' T (z,y))G T (;z')\ 

Since for fixed hj^O, 

u h = i(G T (-,y) - G' T (;y)) + (G' T (z,y) - G' T (z',y))G T (;z) 

+ (G' T (z',y)-G' T (z,y))G T (;z') 

does inherit the integrability properties of Gt( - ,2/) and G' T (-,y) from Corol- 
lary 2.2, we now may conclude that Uh £ L 2 (Z d ), and therefore Uh = 0, that 
is, 

l(G T (x,y) - G' T (x,y)) + (G' T (z,y) - G' T (z' ,y))G T {x,z) 

+ (G' T (z',y)-G' T (z,y))G T (x,z') = 

for every x £ Z d . Since by Lemma 6, G^(x,y;-) is continuous in o(e), we 
learn that Gx{x,y; •) is continuously differentiable w.r.t. o(e) and that (2.15) 
holds. 
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We set for abbreviation 

G T (x,z) - G T (x,z'), 
G T (z,y)-G T (z',y), 

G T (z, z) + G T (z', z') - G T (z, z') - G T (z',z) 



G T (x,e) 
(5.10) Gr(e,y) 
G T (e,e) 

Step 2. Proof of 







(5.11) 



da(e 

d 
da(e 



■G T (x,e) = -G T (e,e)G T (x,e), 
■Gr(e,y) = -G T (e,e)G T (e,y). 



This is a consequence of (2.15) for y = z, z': 




da(e 



■(G T (x,z)-G T (x,z')) 







-G T (x,z) 







-G T (x,z) 



da(e) da(e) 

( = 5) -(G T (x, z) - G T (x, z'))(G T (z, z) - G T (z>, z)) 
+ (G T (x, z) - G T (x, z'))(G T (z, z') - G T (z', z')) 
= -(G T (z, z) + G T (z', z') - G T (z, z') - G T (z, z')) 
x (G T (x,z)-G T (x,z')) 

and for x = z, z', respectively. 
Step 3. Conclusion. 
Note that Corollary 2.3 implies 

(5.12) |G T (e,e)|<l. 

The combination of (5.11) with (5.12) yields 







da(e 



■G T (x,e) 



\G T (x,e) 







da(e 



■G T (e,y) 



\G T (e,y)\. 



Since a(e) is bounded, this also yields 
sup|G T (x,e)| ~ \G T (x,e)\, 

a(e) 



sup|G T (e,y)| ~ \G T {e,y)\, 

a(e) 



which is nothing but (2.16). 
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5.3. Proof of Lemma 2.4- We recall that e = [z, z'], z' = z + e.;. 
Step 1. Proof of (2.12). 

We first give a heuristic argument for (2.12) based on the defining equation 

(5.13) T~ Vt(x) - (V* • A(V0 T (x) + 0)0*0 = 0. 

Differentiating (5.13) w.r.t. a(e) yields as in Step 1 of the proof of Lemma 
2.5 

- 1 d ^ T -(x)-(v*.AVp?-)( X ) 



da(e) \ da(e) 

(5.14) 

- {Vi<fr{z) + £i)(6(x -z)- 5{x - z')) = 0. 
Provided we have S L 2 (Z d ), this yields by definition of G T 

'">) = -(V^t(z) + fc)(G r (x,z / ) - Gr(x^)), 



da{e) 
which is (2.12). 

In order to turn the above into a rigorous argument, we need to argue that 
4>t{x) is differentiable w.r.t. o(e) and that J^fy £ L?{7L d \ Starting point is 

the representation formula from Step 2 of the proof of Lemma 2.6, that is, 

(5.15) <fr{x) = [ G T (x, y)V* ■ (A(y)O dy. 

Jz d 

Combined with Corollary 2.2, (5.15) and (2.15) in Lemma 2.5 show that 
4>t{x) is differentiable w.r.t. a(e). We may now switch the order of the 
differentiation and the sum as follows: 

r (x) = -V Zi G T (x,z)t i 



da(e) 

(5.16) - I V Zi G T (x,z)V Zi G T (z,y)V* ■ (A(y)Ody 

Jz d 

= - y 2i G T {x,z) (^ + J y Zt GT(z,y)y*-(A(y)Q dy) , 



since Gt{-,z) E L 2 (Z d ) by definition of the Green's function, Gt{z,-) £ 

d<f>T 
da(e) 



L x (Z d ) by Corollary 2.2 and A is bounded. This proves that G L 2 (Z d ). 



Step 2. Proof of 

(5.17) sup|^r(x)| < \<fr(x)\ + (\ViMz)\ + l)\V Zi G T (z,x)\, 

o(e) 



(5.18) sup 

o(e) 



d(f> T {x) 



da(e) 



<{\Vi<t>T{z)\ + l)\V Zi G T {z,x) 
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We argue that it is enough to prove (2.14). Indeed, the combination of (2.12), 
(2.16) and (2.14) with the boundedness of a implies (5.17) and (5.18). In 
order to prove (2.14), we proceed as follows 

d4> T \ , , d(f)T , v 84>t , ,s 

Vl d^)) {z) = Me) {z) ~Me) {z) 

{2 = ] {ViM*)+Si)(GT(z,z) - G T (z,z')) 

- (ViMz)+Zi)(GT(z',z)-G T (z',z')) 

(5.19) 

= (Vi<M*) + &) 

x (G T (z, z) - G T (z, z') - G T (z', z) + G T (z',z')) 
= (ViMz) + &GT(e,e), 
where we used the abbreviation 

G T (e,e)=G T (z,z) - G T (z,z') - G T (z' ,z) + G T (z' ,z'). 
Recalling that Corollary 2.3 implies 

G T (e,e)<l, 

inequality (2.14) follows now from (5.19) and the boundedness of a. 
Step 3. Proof of (2.13). 
For n > 0, the chain rule yields 

8a(e) 1 + M J 8a(e) ' 

Using (5.17) and (5.18), this implies 

d[4) T {x) n+l 



sup 

a(e) 



dale) 



(\Mv)\ + (\ViMz)\ + W Zi G T (z, x )\T 



x((|V^t(^| + 1)|V Zi G t (z,x)|), 
which turns into (2.13) using Young's inequality. 

5.4. Proof of Lemma 2.6. We first prove the claim for Gt and deduce 
the result for <fir appealing to an integral representation using the Green's 
function. 

Step 1. Properties of Gt- 

The product topology is the topology of componentwise convergence. 
Hence, we consider an arbitrary sequence {a^}^^ C A a p of coefficients such 
that 

(5.20) lim a u (e) = o(e) for all edges e. 

isfoa 
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Fix y G Z, d ; by the uniform bounds on GT(-,y,a u ) from Corollary 2.2, we 
can select a subsequence v' such that 

(5.21) ut(x) := lim Gt(x,?/;cv) exists for all x G Z d . 

i/'too 

It remains to argue that ut(x) = GT(x,y;a). Because of (5.20) and (5.21), 
we can pass to the limit in (T~ 1 Gt(-, y; a u >) + L a ^, 6?t(-, y; a v >)){x) = 5(x — y) 
to obtain 

(5.22) (T _1 it T + LatiT)(x) = 5(z - y) for all x G Z d . 

Moreover, the uniform decay of Gr(-,y;a v ) from Corollary 2.2 is preserved 
in the limit, so that ut G L 1 (Z d ) C L 2 (Z d ). Note that Riesz's representation 
theorem on L 2 (Z rf ) yields uniqueness for the solution of (5.22) in L 2 (Z d ). 
Hence, we conclude as desired that ut(-) = G*T(-,y;a). Borel measurability 
of GT(x,y;-) in the sense of Lemma 2.3 follows from continuity w.r.t. the 
product topology, cf. [10], Satz 14.8. 
Step 2. Properties of 4>t- 

Corollary 2.2 ensures that G T (x, •) G L x {1 d ) for all x G Z d and one may 
then define a function (pT by 

(5.23) Mx)= I G T (x,y)V*-(A(y)Zi)dy. 

Since Gt(- + z, - + z) has the same law as Gt(-, •) by uniqueness of the Green's 
function and joint stationarity of the coefficient A, c/>t(' + z) has the same 
law as (pT- This shows that 4>t is stationary. In addition, <J)t is a solution 
of (2.3) by construction. Hence, by the uniqueness of stationary solutions 
of (2.3), (j)T = 4>T almost surely, so that by the measurability properties we 
may assume 4>t = 4>T- 
Introducing for R > 1 

<I>t,r{x) ■= I G T (x, y)V* • (A(y)d) dy, 
J\v\<R 

one may rewrite (5.23) as 

(5.24) Mx)= lim <Pt,r{x). 

R-^-oo 

From Step 1, 4>t,r{x) is a continuous function of a since GT{x,y) is and the 
formula for 4>t,r(x) involves only a finite number of operations. Note that 
Corollary 2.2 implies that 

lim sup / GT(x,y;a) dy = 0. 

Rtoo adA aP J\y\>R 

Hence, the convergence in (5.24) is uniform in a and the continuity of <Pt,r in 
a is preserved at the limit. Therefore, 4>t (and continuous functions thereof) 
are continuous with respect to the product topology, and hence Borel mea- 
surable. 
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5.5. Proof of Lemma 2.7. We first sketch the proof in the continuous 
case, that is, with Z rf replaced by M. d . 
Step 1. Continuous version. 

Starting point is the defining equation (2.3) of the corrector 4>t in its 
continuous version, that is, 

(5.25) T~ Vt - V • A(V(j)T + = in R d . 

We multiply (5.25) with 4>^ +1 and obtain by Leibniz' rule: 
= T~ Vt +2 + (-V • AiVfa + 0)<Pt +1 
= t- 1 ^ 2 - V • (^ +1 A(V<Ar + 0) + V( At +1 • A (V<h + 

(5.26) 

= T- 1 ^ +2 -V-(^ +1 A(V<Ar + e)) 

We then take the expected value. Since the random fields A and 0r are 
jointly stationary, and thus also <^ +1 ^4(V 4>t + £), we obtain 

{T~ l r T +2 ) + (n + lXffiVfa ■ A(V<p T + 0) = 0, 

and therefore 

(ffiV(fr-A(V(f>T + £))<0 

since n + 2 is even. By the uniform ellipticity of A and since (f> j. > (n is 
even) and |£| =1, this yields the estimate 

<0?Wt| 2 > < (#|V0 T |). 

Applying Cauchy-Schwarz' inequality in probability on the right-hand side 
of this inequality yields the continuum version of (2.17), that is, 

(^|V0 T | 2 )<(^). 

We now turn to our discrete case. 
Step 2. Discrete version. 

We need a discrete version of the Leibniz rule V • (fg) = /V • g + V/ • g 
used in (5.26). Let / G Lf oc (Z d ) and g G Lf oc (Z d ,R d ), then this formula is 
replaced by 

d 

V* • (fg)(z) = ^(f(zMz)]j - f(z - ej )[g(z - e,)],) 
i=i 

(5.27) 

cZ 

= /(z)V* • g(z) + V*f(z)[g(z - e,)];. 
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We also need a substitute for the identity V0^ +1 = (n + l)0yV0y used in 
(5.26). This substitute is provided by the two calculus estimates 

(5.28) (4> n+1 - r +1 )(4> - <t>) > (4> n + <P n )(4> - <P) 2 , 

(5.29) |^+i_^n+i|<(^ + ^»)|^_4 

For the convenience of the reader, we sketch their proof: by the well-known 
formula for cj) n+1 — <ft n+1 , they are equivalent to 



m=0 

By homogeneity, we may assume (fi = l, so that the above turns into 



Yl m ~ 1 + 



m=0 

The upper estimate is obvious by Holder's inequality since n is even. Also for 
the lower bound, we use the evenness of n to rearrange the sum as follows: 

n i i 1 

E ( ^ m =2 1+ 2 (1 + 20 + (/,2) + 2 (/,2(1 + 2</, + </,2) + '" 

m=0 

+ ^-2 (1 + 2 ^ + ^2 )+ l n 

After these motivations and preparations, we turn to the proof of Lemma 
2.7 proper. With f(z) :=<^; +1 (z) and g{z) := A(V<fr + ( 5 - 27 ) turns 

into 

v*-(^ +1 (^(V0 T + O(^)) 

d 

= ^ +1 (z)V*-A(V<At + 0(^) + E V ^t +1 (^) ^(V0T+e)(^-e,-)b • 

7=1 ' ' 

=a(z-e:j,z)(Vj4>T(z-ej)+$j) 
=a{z-e j ,z)(V*<t> T {z)+iii) 

Hence, 

-^ +1 (z)V*-A(V^ r + C)(^ 
d 

(5.30) = ^V*^ +1 (z)a(z - ei ,z)(V^ T (z) +&) 

j'=i 

-v*-(^ +1 (^(V0t+6(^))- 
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Multiplying (2.3) with 4>t +1 (z) and using (5.30) emulate (5.26) and yield 
= T- l ^+\z) - V* • (<^+\z)A(y<h + £)(*)) 

(5.31) 

d 

+ E V*^ +1 (z)a(z - ej ,*)(V;<Mz) + 

Taking the expectation of (5.31) and noting that <^. +2 > 0, we obtain as for 
the continuous case 

(Ta(z - ej ,z)V*^ +1 (z)V*0 T (z)^) 

(5.32) 3 

<^a(z- ei ,z)|V*^ +1 (z)|y 
On the one hand, we have 

d 

^a(z-e„z)V*^ +1 (z)V^ r (z) 
i=i 

d 

(5.33) = ]Ta(z - e i5 z)(^ +1 (z) - ^ +1 (z - e,))^^) - ~ 

i=i 

(5.28) d 

> £>£(*) + 0t(* - *MM*) - Mz - e,)) 2 - 

On the other hand, we observe 

^a(z- ei ,z)|V*^ +1 WI 

d 

(5.34) = E°( z " e i5 2)|<^ +1 W - r T + \z ~ 

i=i 

(5.29) d 

Now (5.32), (5.33) and (5.34) combine to 
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d 

< T,m(z) + M* - e,))\Mz) - <fr{z ~ ej)\). 
i=i 

By stochastic homogeneity, this reduces to 

^E(0T(ei) + <ft(0))(M*,) - <M0)) 2 ^ 

An application of Cauchy-Schwarz inequality yields 

^(^)+^(o))(M^)-Mo)) a ^<^(^)+^(o))y 

A last application of stochastic homogeneity gives as desired 

^(o) Y,((M*j) - Mo)) 2 + (Mo) - M-*i)) 2 )] Z <^(o)>- 

5.6. Proof of Lemma 2.10. The proof relies on a doubly dyadic decom- 
position of space. First note that by symmetry, 

/ hr(z)hr(z — x) dz = / hr(z)hT(z — x) dz 

J\z\<\z-x\ J\z\>\z-x\ 

1 f 

> — I hx{z)hT{z — x) dz. 

2 Jz d 

Hence, it is enough to consider 

/ / hT(z)hx{z — x)dzdx. 

J\x\<rJ\z\<\z-x\ 

In the three first steps, we treat the case d > 2. We then sketch the modifi- 
cation for d = 2 in the last step. Let R ~ 1 be such that (2.27) holds with a 
constant independent of R for all R > R/2. 
Step 1. Proof of 

(5.35) / / h T (z)h T (z-x)dzdx<R 2 for R > 2R. 

JR<\x\<2R J\z\<\z-x\ 

Let JV e N be such that R<2~ N R<2R.We then decompose the sum over 
\z\ < \z — x\ into three contributions: R/2 <\z\, a dyadic decomposition for 
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R < \z\ < R/2 and a remainder on \z\ < R. More precisely, 



/ / hx(z)hT(z — x) dzdx 

JR<\x\<2R J\z\<\z-x\ 

= / hT(z)fiT(z — x) dzdx 

Jr<\x\<2R Jr/2<\z\<\z-x\ 



R<\x\<2R J R/2<\z\<\z-x\ 
N 

+ y / / hT{z)fiT{z — x) dzdx 

n=1 JR<\x\<2R J {2-( n + 1 'l R<\z\<2- n R}n{\z\<\z-x\} 



+ / h,T(z)hT{z — x) dz dx 

J R<\x\<2R J {\z\<2-( N + 1 ) R}n{\z\<\z-x\} 

< / Kt{z)}it{z — x) dz dx 

J\x\<2R J Rl2<\z\<\z-x\ 



=h 

+ y / / hr(z)hT(z — x) dzdx 

~[ JR<\x\<2R J2~("+ 1 )R<\z\<2-"R 

l v ^ , 

=h(n) 

+ / fiT(z)hT(z — x) dz dx . 

J R<\x\<2R J \z\<R 

y „ ' 

=h(N) 

We use Young's inequality, a dyadic decomposition of {\z\ > R/2}, and the 
assumption (2.27) to bound I\: 

h<\( ! h T {zfdz+ [ h T {z-x) 2 dz) 

2 \JR/2<\z\ JR/2<\z-x\ J 

00 r (2-27) 00 / 1 \ k 

In order to bound l2(n), we will use the following fact 
(5.36) (|x| > R and \z\ < \R) (\z -x\> \R). 

We have by Cauchy-Schwarz inequality 

J 2 (n) < ( / / h T {zfdzdx 

\J\x\<2R J2-( n + 1 )R<\z\<2-"R, 



/ / ^t(-2 — a;) 2 dzdx J 

/r<|:e|<2.R J|z|<2-™ J R / 
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<(R a h T (zYdz 

1 2-(™+!) R<\z\<2~ n R 



(2.27) 
< (2-"R) 2 - d 



1/2 

z — x) 2 dzdx 




hx(z — x) 2 dxdz 

'\z\<2~ n RJ R/2<\z-x\<5R/2 

1 V 

(2 | 7) / R 2 - d dz = (2~ n R) d R 2 - d 

J|z|<2-™_R 



<2-»/2 2 . 



We proceed the same way to bound Is(N). Recalling that R > 2R ~ 1, it 
holds that \z\ < i?=>- |z| < i?/2. Hence, we are in position to use (5.36) and 
we obtain 



h(N)<( [ [ h T (z) 2 dzdx [ [ h T (z-x) 2 dzdx] 

\J\x\<2RJ\z\<R Jr<\x\<2RJ\z\<R J 

<(W h T (z) 2 dz 

\ J\z\<R 



1/2 

Jit{z — x) 2 dzdx 

\z\<R J R<\x\<2R J \z\<R 

v ' * v ' 

(2.28) (5.36) f f „ 

< i < / / hr(z — x) \ dxdz 

J\z\<R J R/2<\z-x\<5R/2 



(2 f } / R 2 ~ d dz~R 2 ~ d 

J\z\<R 



<R. 



Since £^° =1 2~ n R 2 ~ R 2 and |{|x| < 2 J R}|i? 2 - d ~ i? 2 , the bounds on h, I 2 (n) 
and Is(N) imply the claim (5.35). 
Step 2. Proof of 

(5.37) / / h T {z)h T (z-x)dzdx<\. 

J\x\<m J\z\<\z-x\ 

This time, we decompose the sum over \z\ < \z — x\ in two contributions 
only: \z\ < R and R< \z\. We then obtain 

/ / hr(z)hr(z — x) dzdx 

J\x\<4R,J\z\<\z-x\ 
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/ hr(z)hT(z — x) dz dx 

\x\<4R J R<\z\<\z-x\ 



-I'l 



+ / hT{z)hr(z — x) dz dx . 

'|x|<4R J {\z\<R}n{\z\<\z-x\} 



Proceeding as for I\ in Step 1 using (2.27) yields 

For I' 2 , we use Cauchy-Schwarz inequality, (2.28), and i?~ 1: 
T 2 <( [ [ h 2 T (z)dzdx] 1 ( [ [ h 2 T {z')dz' dx\ ' < 1. 

\J\x\<iR J\z\<R J \J\x\<4R J\z'\<5R / 

This proves (5.37). 

Step 3. Proof of (2.30). 

It only remains to use a dyadic decomposition of the ball of radius R into 
the ball of radius R and annuli of the form 2~ k R <\z\< 2~ k+1 R, as follows. 
Taking M such that 2R < 2~ M R < 4R, it holds that 



hT(z)hT(z — x) dzdx 

\x\<RJ\z\<\z-x\ 



hT(z)h,T(z — x) dzdx 

\x\<2~ m RJ\z\<\z-x\ 



(5.37) 
< 1 



M 

+ / / / hr(z)hT(z — x) dzdx 

^—J J 2 ™-M-iR < \x\<2 n - M R J\z\<\z-x\ 



(5.35) 

< (2"- M i?) 2 



M 



<i + R 2 Y,^ n ~R 2 , 



which proves (2.30). 
Step 4. Proof of (2.29). 

For the case d = 2, we use the same strategy as for d > 2. The bounds on 
l2{n) and Is(N) are the same as for d > 2. However the estimate for I\ is 
slightly worse. Indeed, we split the dyadic sums 2 k R< \z\ < 2 k+1 R into two 
categories in order to take advantage of the fast decay in (2.26): the first class 
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is for k such that 2 k R< y/T and the other class for k such that 2 k R> \/T. 
More precisely, setting 1(R, T) := {k G N : 2 k ~ l R < VT}, we have 

hT{z)fiT(z — x) dz 

R/2<\z\<\z-x\ 
Young f 

< \ h T {z) 2 dz 

JR/2<\z\ 



= V / h T (z) 2 dz 

k= _ 1 J2 k R<\z\<2 k + 1 R 

= h\(z)dz+ y~] / h 2 T {z)dz 

kex(R,T) j2k - 1R <^ 2kR " keN\x{R,T) j2k - lR <^ 2kR 



(2.26) (2.26) 

< max{0,ln(v / T^ 1 )} < £ 2- 2fc <l 

fcew 

< max{l, ln( VTiT 1 )}, 
which gives the extra factor in (2.29). 

APPENDIX: HEURISTICS FOR (1.13) AND (1.14) 

Let (pi and 4>t,i denote for i G {1, . . . ,d} the solutions of (1.11) and (1.12), 
respectively, with £ replaced by the ith unit vector of M. d . We claim that 

d d 

(A.l) J2J2 vas [J2^- (A-(A))ei + 2 ej - V^i] = dvar[a]^£, 
i=i j=i 

d 

(A.2) V^ r ,i - V^| 2 ) = var[a]T- 2 G 2 ., 

i=l 

where Gt denotes the fundamental solution of the constant coefficient opera- 
tor — A. We also denote by G the fundamental solution of the Laplacian. 
Since 

f T 2-d/2^ ford<4, 
^ G r~ I hiT, for d = 4, 

[ 1, for d > 4, 

and 

(1.13) and (1.14) follow from (A.l) and (A.2), that we prove now. 
Step 1. Argument for (A.2). 
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Since 

(A.3) - A$r - $ = -T" 1 ^, 

one has 

(A.4) (|V(0 T - 4>)\ 2 ) = -T-\4> T {h - 4>)). 

Rewriting (A.3) in the form 

T~\4> T -$)- A^t -4>) = -T- 1 ^ 

yields the formula 

(A.5) {4> T - 0)(O) = -T- 1 J3Gr(x)^(x). 

X 

Using (A.5), (A.4) turns into 

(A.6) (|V(0 T - 0)| 2 > = -T- 2 ^G T (x)(^ T (O)0(x)). 

X 

Expressing now </>Ti(0) and 4>i(x) in terms the Green's functions 2 Gt and 
G, 

^T, i (0) = ^G T (x , )V*-(^(x / )e i ) 

x> 

= -Y,V i G T (x')(a i (x')-(a)), 

x> 

4> i {x) = Y J G(x-x')V* ■(A(x')ei) 

x> 

= -J2^iG(x-x")(a i (x")-(a)), 

x" 

and using the independence of di(x') and cii(x") for x' ^ x" , we get 
(4r,i(P)Ux)) = ^V i G T (x')V i G(x-x / )((a i (x / ) - (a)) 2 ). 

x' 

Hence, 

d 

5^(0r,i(O)0i(a:)> = var[a] ^ VG^') • VG(x - x') 



i=l 



■ v&r[a]GT(x) 



Attention should be paid here to turn this into a rigorous argument since G is not in 
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since -AG(x) = 5(x). Combined with (A. 6), this proves (A. 2). 
Step 2. Argument for (A.l). 
Using the Green's function, one has 

Mx) = Y^G{x- x')V* ■ ((A - (A))e t )(x') 

x> 

= -J2ViG(x-x')(ai(x')-(a)), 

x> 

and therefore 

Vfaix) = -^WiG(x - x'){ ai {x') - (a)). 

x' 

Hence, denoting by Aij the argument of the variance in (A.l), one has 
A? := X^ e -? ' ei(ai(i) - (a)) + 2ej ■ Vcpi(x))r] L (x) 

X 

= Y,^2( a i( x> ) - ( a )) e j • ( S ( x - x')ei - 2VViG{x - x')) VL (x). 

x x' 

Using the independence of the aj, one obtains for the variance 

var[Aj] = var[a] ^ ^ ^ e r (5(x - x')ei - 2VViG(x - x')) 

x x' x" 

x ej • (8(x" - x')ei - 2VV l G{x" - x')) 
x r] L (x)r] L (x"). 

Rearranging the terms yields 

var[Aj] = var[a] ^^<5(j - i)(S(x - x) - 4ViViG(x - x'))t] L (x)t] L (x / ) 

X X 1 

+ var[a] ^ ^ ^} L {x)r] L {x") V j V,G(x - x') Vj V;G(x" - x') . 

X x" x' 

" „ " ' 

=-V I V l E x ' G{x-x')V j V j G(x"-x>) 
Summing in j and using that — AG(x) = S(x), this turns into 
d 

^var[Aj] = var[a] ^^(<5(x - x') - 4ViViG{x - x'))i] L (x)r] L (x') 

j=l x x ' 

+ var[a] 4r ?L (x)r ?L (x ,/ )V i V i G'(x - x") 

X x" 

= var[a] 5(x — x')r]L(x)r]i(x') 

x x' 
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= vai[a]^2f] L {x) 2 , 

X 

from which we deduce (A.l). 
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